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NOMENCLATURE 


National Meteorological Center octagonal-grid coordinate systen. 
rectilinear coordinate system for SPUOK, SCAWUEN, RAIN, SHASTA, 
and SHASTA (km). 

rectilinear coordinate system for TEMPEST (km). 

altitude (km). 

time-step index. 

(1519552555) = phase-space point in SPUOR at time step j. 
velocity recursion-relaticnship correlation parameter between 
time steps j and i. 

distance (km). 

total fission-equivalent yield of a SPOUR source cloud (kt). 
time at step j measured from the Spru0R particle's birth (h). 
time-mesh width (h). 

total simulation time in a SPOUR calculation; or cbservation time; 
or maximum threat-time (h). 

total number of time steps in a SPOUK-particle trajectory; or 
total population with lwbD > D; or overall radial-binning param- 
eter for the effects-tallying grid in SCAWUEN. 

total number of source particles in a SPUOK calculation. 

average debris-cloud centroid velocity (km/h). 

orthogonal wind components in the NCAn wind data oase (km/h). 
rain area fraction. 

average rain area fraction. 

precipitation-cell area distribution parameters for TEMPEST. 
maximum precipitation-cell area (km*). 
semimajor axis of a precipitation cell (km). 
semiminor axis of a precipitation cell (km). 
maximum value of a during cell's lifetime (km). 
total lifetime of a cell (h). 
precipitation-cell aspect ratio = b/a. 


location of a precipitation-cell center in TEMPEST. 


V = (Vy,Vy) = precipitation-cell velocity in TEMPESi (km/h). 
¢ = angle between Y-axis and v. 
6 = angle between cell minor axis and v. 
port dé. 
a, = (X51V 522703) = phase-space point in SCAWUEN. 


(Xo,¥q) = location of a precipitation-cell center in SCAWUEN. 
(X5 57125): = location of a SPOUR particle within a cell in SCAWUEN. 
R. = maximum rain intensity at the center of a precipitation cell (mm/h). 
Voy Vy = within-cell intensity distribution parameters. 
0, = Yya. 
Op = Yyb. 
a(Z) = vertical-intensity distribution. 
Kk = local precipitation rate (mm/h). 
dX = seavenging coefficient (h71), 
We = total debris-particle weight scavenged during time At from a given 


SPOOR particle at :. 


J 
Ww, = unscavenged weight of a SPOUR particle at a. 
v, = viz) = horizontal transport velocity at altitude z for a scavenged 


debris particle (km/h). 
tp = time for a scavenged particle to reach the ground (h). 
AW = ground-deposited weight during time step j from SPUOK particle m. 


jm 
CX imeY jm) = ground location of AW jm 
(Xjos¥ jo) = location of the centroid of ground-deposited particles during 
time step j. 
P jm = radial distance from the centroid of ground-deposited particle m 


at time step j (km). 

i ¥, = first moment about the centroid of the radial distribution of 
| ground-deposited particles at time step j (km). 

of = second moment about the centroid of the radial distribution of 
ground-deposited particles at time step j (km). 

| ; kK = number of angular divisions for the effects-tallying grid in 


SCAWUEN; or the amplitude of the scavenging-coefficient function. 


Fs outer radial-binning parameter for the effects-tallying grid in 
SCAWUEN. . 
‘Ty = radius of the nth annular ring in the effects-tallying grid 


in SCAWUEN (km). 


vi . ; 


largest of the P jm (km). 

area of tally-bin k at time step j (km). 

deposited-debris intensity for a tally-bin within which point k 
is located at time step j (particles/km“). 

free-field infinite whole-body dose (rem). 

dose rate at time t after detonation from deposited debris (rem/h). 
time after detonation thut deposition occurs (nh). 

IwkED at point k from deposition during time step j (rem). 

total Iwbb at point k for all time steps (rem). 

average precipitation rate for observation time T (mm/h). 
probability that RAF > HAF, given that RAF # 0. 

probability that observable precipitation occurs during observa- 
tion time t. 

probability that Kar # 0 during observation time T and that RAF 
> KAR. 

occurrence probability. 

maximum population with 1lwWbD > Dy. 


naximum threat-distance (km). 


vil 


HALNUUT ASSESSMENT: THE ACKA SYSTEM AND SUMMARIES OF STIMULATION RESULTS 


by 
Clayton W. Watson, Sumner barr, anc Ray &. Allenson 


AbSTRACT 

We developed a generalized, three-dimensional, integrated computer code 
system to estimate collateral-damage threats from precipitation-scavenging 
("rainout") of airborne deoris~clouds from defensive tactical nuclear 
engagements. This code system, called ACKA for Atmospheric-Contaminant Rainout 
Assessment, is based on Monte Carlo statistical simulation methods that allow 
realistic, unbiased simulations of probabilistic storm, wind, and precipitation 
fields that determine actual magnitudes and probabilities of rainout threats. 
Detailed models (or data bases) are included for synoptic-scale storm and wind 
fields; debris transport and dispersal (with the roles of complex flow fields, 
time-dependent diffusion, and multidimensional shear effects accounted for 
automatically); microscopic debris-precipitation interactions and scavenging 
probabilities; air-to-ground debris transport; local demographic features, for 
assessing actual threats to populations; and nonlinear effects accumulations 
from multishot scenarios. we simulated several hundred repres2ntative shots 
for West European scenarios and climates to study single~shot and multishot 
sensitivities of rainout effects to variations in pertinent physical variables. 
One interesting result was the consistent behavior of rainout effects versus 
scavenging efficiency A, rain rate R, and shot yield Y: effects dropped off 
rapidly at A < 0.1 h71, R ¢ 0.1 mm/h, or ¥Y S 0.1-0.4 kt. More important, we 
demonstrated broad, shallow, maximum-threat regimes for } 2 0.5 h-1, 0.1 . K < 
1 mm/h, and 0.4 s 1's 5 kt. Other inherent features of rainout determined for 
a spectrum of scenarios were maximum threat-times and distances, occurrence 
probabilities in the west German climate, and comparisons between rainout 


effects and fallout or prompt collateral-damage effects. 


I. INTRODUCTION 

Defensive tactical nuclear engagements in adverse weather involve the 
possibility that precipitation scavenging ("rainout"*) of airborne nuclear 
weapon debris could deposit highly radioactive material on the ground. If such 
deposition ocecurs before dispersal in the atmosphere can dilute the debris 
cloud, this "wet fallout" could pose a delayed collateral-damage threat to 
downwind populations. 

Depending on their likelihood and magnitude, these threats may have 
important implications for future US nuclear-weapon design and employment 
doctrine. Rainout effects, for example, could be of particular significance 
when weapons of <10-kt fission yield are used in an airburst mode. In this 
mode, a device is detonated at an altitude that precludes direct interaction of 
the fireball with the ground surface, and littie or no surface material mixes 
with the weapon vebris. In this case, the radioactive debris is attached to 
aerosol-like particles so small that they remain aloft, and conventional, 
gravity-driven, “dry fallout is essentially eliminated. Thus, for a weapon 
that is airburst to maximize useful military effects and eliminate dry fallout, 
rainout becomes the only remaining source of delayed collateral damage. The 
question is then whether, or under what circumstances, rainout effects might 
negate the collateral-damage advantage of airbursts and how the problem, if it 
exists, can be minimized. 

The Los Alamos Scientific Laboratory (LASL) has an ongoing study program 
to assess rainout effects and implications. lhe have developed a generalized, 
three-dimensional, integrated, computer code system for this purpose and have 
used this system to estimate rainout effects for several hundred representative 
nuclear shots in west European scenarios and climates. In this report, we 
describe our rainout assessment system and its underlying physical bases, and 
present numerical results. Our calculations have given us considerable insight 


into the character and magnitude of rainout threats. 


II. BACKGROUND 
Rainout has posed a persistent and controversial problem in the US 


*In this report, we will use the term "rainout" both in a generic sense 

to imply any and all precipitation-scavenging events and in the more 
conventional way to specifically denote in-cloud scavenging events. The work 
reported here deals primarily with this latter phenomenon, which is generally 
recognized as the “.minant source of potential collateral-damage threats.'- 


nuclear~weapon community for several years. Historically. some confusion seems 
to have developed as to what the relative roles and importance of the many 
physical elements involved are and what we should calculate. and why. 

we believe that uncertainties as to what the roles of the driving elements 
are can be traced to overall modeling approaches that have been too narrow in 
outlook and to the resulting, often subtle, rigidities built into the 
assessment codes. it is reasonable, for example, to make conservative 
estimates for the elements of a problem and then combine them into upper-limit 
estimates of effects, provided the measures of conservatism are adequately 
understood and ihe results are truly bounding; then, if the results are within 
acceptable limits, the problem goes away. Un the other hand, if such estimates 
show that a potential problem stilt exists (as with rainout), the assessment 
procedures and tools must be upgraded until uncertainties diminish enough for 
decision makers to understand and use tne model results with confidence. 

fainout assessment is an especially insidious task. fundamental 
complexities in the atmospheric transport of debris clouds, in the driving 
meteorological phenomena, and in the precipitation processes themselves make 
the analytical problems difficult, involving space and time scales that span 
many orders of magnitude. moreover, rainout threat levels and damage 
probabilities are determined by intricate interactions and tradeoffs among 
several physical, logistic, and demographic factors; the problem is tnus 
intractable in broad terms. Intuition is also often misleading, and 
generalized ccnciusions usually are too easily begged somewhere in the 
analytical process. 

we recognized in mid-1473 that rainout assessment had reached an impasse 
because of the avoove difficulties and because we lacked a sufficiently general 
and flexible modeling framework for systematically examining the physical 
elements of the problem to place them in proper perspective. A totally new 
approach was needed that emphasized pragmatic goals while still recognizing and 
effectively medeling tne relevant physical phenomena. The approach also had to 
be directed explicitly to acquiring sensitivity information needed to estimate 
actual rainout effects for realistic scenarios, and it had to establish 
confidence bounds within which these estimates could be used. 

A major goal of our work has been to develop such an approach and the 


analytical tools to implement it. ‘the basis of our scheme is an overall 


modeling framework designed to be general and flexible enough to accommodate 
continual upgrading of model components to facilitate a progression through a 
hierarchy of increasingly sophisticated analyses, as required, to achieve 
acceptable accuracy in the results. This approach depends on continual 
information feedback between model development and technical analysis (Fig. 1) 
te identify and evaluate relevant factors, avoid irrelevant or secondary 


issues, and avoid the “intuition” and "generalization" pitfalls nentioned 


earlier. 


IIl1. PROBABILISTIC NATURE OF RAINOUT 

Rainout is inherently stochastic--a direct result of the capricious nature 
of local wind, storm, and precipitation patterns. The most important 
stochastic variables relate to the time and space structure of precipitation 
celis and to turbulent air motions that determine diffusion and meander of 
aebris clouds and entrainment of debris into precipitating systems. This 
character of rainout means that rainout collateral-damage effects are 
fundamentally unpredictable except in probabilistic terms. Upper-limit 
estimates thus do not provide adequate insight into rainout threats. Because 
rainout damage is a relatively unlikely event, we are interested in unlikely 
combinations of the relevant individual physical events. Practical cases of 
interest lie in the "tails" of the governing probability density functions 
(PDFs), and this fact, plus the fact that the analytical problems are complex, 
forces us so far out on the PDF tails in attempts to be "conservative" that the 
resultant overall rainout situations are hardly credible. In other words, to 
impose extreme or bounding conditions on these probabilistic variables without 
begging the question becomes impossible, particularly since it is often 


difficult even to know what constitutes a conservative assumption. 
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Fig. 1. ; 
Dynamic information feedback. 


Unfortunately, "average" or "typical" estimates of rainout effects are not 
very informative either, because rainout is itself atypical.* Real 
rainout-threat situations are never average, nor are any resulting rainout 
effects, and the less likely events are of greater interest. For example, the 
consequences of a low-probability rainout event that happens to occur over a 
large city can be more significant than a high-probability event where no 
Significant population centers lie downwind from the shot. Ina given 
instance, furthermore, the total hazard represented by the combined 
probabilities of several unlikely events may pose a problem even though each 
individual threat is innocuous. 

The key point here is that not only must we estimate hazards and 
consequences for representative specific cases, but, more important, we must 
also estimate the associated probabilities. In this process, it. is also 
essential that we be mathematically correct in combining the probabilities of 
the pertinent individual elements so that we obtain meaningful overall rainout 
assessments. For these reasons, we have based our analytical approach on Monte 
Carlo statistical simulation techniques that allow us to combine the physical 
elements of the problem relatively easily in a statistically unbiased, 
multidimensional fashion; because individual elements are combined 
probabilistically, essential stochastic variabilities are accounted for 
naturally. 

The logical core of our simulation system is a three-dimensional, Monte 
Carlo, debris-transport model; other elements are modeled in a similar fashion 
based on integral inputs in the form of PDFs extracted from the literature and 
from available data bases. Collateral-damage estimates for specific 
battlefield scenarios depend on the collocation of debris, storm systems, 


precipitation, and populations and are available through straightforward, 


linear, statistical tallies. 


IV. REQUIREMENTS FOh A RAINOUT COLLATERAL-DAMAGE EVENT 

The stochastic features of rainout mean that several events must occur at 
the same time and place before significant rainout collateral-damage effects 
can result. These events, along with the principal physical factors that 


determine their distributions in time and space, are listed below. 


#The most frequently occurring rainout-threat level is zero. 


@ bebris must be concentrated enough to pose a threat. Generally, for 


the low-yield devices pertinent to rainout, this condition is met only for a 
few hours and within several tens of kilometers downwind from the shot. 
Griving factors: scenario (particularly yields and number of shots), local 
winds, and atmospheric diffusion rates. 

@ A precipitating storm system must be in the area. kven in the 
relatively wet climate of Western Lurope, most days throughout a year are free 
of large-scale storm sstems that can offer substantial rainout possibilities 


in a given area. briving factors: local climate and storm-system synoptic 


characteristics. 

e Precipitation cells must exist in the area and must be intercepted by 
the debris. kven in very "rainy" parts of a storm, the precipitation patterns 
are cellular in structure; and a specific debris cloud will encounter important 
precipitation cells relatively infrequently. This point requires special 
attention because previous work in rainout has essentially overlooked the 
variability of precipitation on scales comparable to the dimensions of the 
debris clouds and to the separation distances between population centers. This 
precipitation variability is important in determining actual rainout effects. 
Criving factors: mesoscale precipitation morphology, debris and 
precipitation-cell interaction dynamics, and local winds. 

@ Ground-deposition rates must be substantial. both the precipitation 
rate and the scavenging probability per unit precipitation rate must exceed 
certain levels to produce ground doses above a damage threshold. Considerable 
controversy exists regarding scavenging mechanisms and probabilities pervcinent 
to airburst debris. Indeed, there is disagreement as to whether real 
scavenging rates permit rainout collateral damage to occur at all. ven if the 
most "dangerous" scavenging mechanisms dominate, however, that in itself does 
not necessarily imply large rainout casualties. Driving factors: scavenging 
mechanisms, debris particle sizes, debris physical-chemical characteristics, 
and precipitation-intensity distributions. 

@ Debris deposition must be over people. We disagree on this point with 
other rainout investigators who use rainout "area" as a measure of rainout 
effects. The concept of area as a measure of collateral damage from delayed 
radiation sources transported through the atmosphere is useful only in certain 


..Special contexts. For example, it can be useful for conventional faliout 


because the "deposition" function (gravity) is constant and we ean, at least in 
principle, predict where a given dose-rate contour will lie. It is then 
reasonable to introduce separately a "damage" function (distribution of people) 
to estimate effects. No such direct relationship is possible for rainout, 
nowever, because both the deposition function (precipitation cells) and damage 
function (people) are distributed discretely and heterogeneously. 


One must therefore distinguish clearly and carefully among concepts such 


as (1) people at risk (because the debris cloud passes overhead and could 
deposit material), (2) "average" number of people under "average" rainout 
depositicn contours, and (3) expected-value probabilities for numbers of people 
exposed to different hazard levels in representative military and 
meteorological situations. The first two items can be compared with prompt- 
and dry-fallout collateral-damage areas, but they are not measures of expected 
rainout damage for a specific situation as is the third item;* nor are they 
related in a useful way to actual rainout effects (except, of course, for the 
special case where the rainout-threatened area is so small that quantitative 
estimates of actual casualty expectations are unnecessary). Driving factor: 
local demugraphy. 

The combined probabilities for coincidence in time and space of the above 
five events (debris, storm system, precipitation, scavenging, and population) 
determine rainout hazards. This combination is a complex mathematical problem 
that cannot be handled by intuition or by multiplying averages together. 
However, when the probability of any one of these events is small for a given 
scenario or if a combination of any subset of these probabilities is small, the 
corresponding overall probability of rainout hazard is small. This point is 
important in gaining perspective on the seriousness of the rainout problem for 


specific contexts. 


V. PHYSICAL ELEMENTS THAT DETERMINE RAINOUT EFFECTS 

In this section, we identify the principal physical elements that 
determine rainout hazards and probabilities, doing so more or less 
chronologically in terms of the debris-cloud history from formation through 


transport, scavenging, and deposition processes to assessment of on-the-ground 


effects. 


*In contrast to gravity-driven fallout, for example, the usual rainout hazard 
is zero. 


A. Scenario. Included here are those parameters that define the configuration 
of the debris clouds at stabilization time (generally a few minutes after shot 
time). These parameters establish the initial conditions or "source" for the 
rest of the problem. 

-The most important factors are number of shots, shot yields, number of 
fissions, and shot locations with respect to each other and to downwind 
populations. Yield directly affects the dimensions of the debris cloud, the 
intensity of radioactivity that may later be rained out, and the vertical 
location of the stabilized debris in the atmosphere. Debris altitude is 
important because the pertinent wind fields and precipitation probabilities 
vary with altitude. (Because precipitation systems are limited to altitudes 
below 10 km, our rainout interest is restricted to devices with yields less 
than a few tens of kilotons.) Number of shots, shot locations, and total 
fissions determine total dose levels and the number of. people threatened at 
various dose levels. 

5. Debris Transport and Dispersal. Debris transport and dispersal establish 
where individual debris clouds will go, how they will grow, and how multiple 


clouds will separate or converge. These are key factors in rainout-effects 


assessment. Debris dispersal (incoherent transport) limits the lifetime of 
rainout threats to one to a few hovrs for a given shot, depending on the yield 
and the damage threshold zonsidered, and determines the width of the threatened 
area on the ground. Debris transport (coherent transport) determines how far 
the cloud will travel during its threat lifetime and, hence, the downwind 
length of the threat area. 

(we emphasize again, however, that the size of this threat area and the 
number of people within it are not measures of collateral damage; they only 
establish bounds. People within this contour are people at risk, i.e., they 
are potential candidates for a rainout event because of the cloud's passing 
overhead, not casualties. Rainout effects can occur only where precipitation 
cells (1) intercept the debris cloud and (2) deposit the debris over people. 
The number of people at risk will usually be large, whereas the number of 
actual casualties will usually be zero and will never be more than a small 
fraction of the people at risk because of the heterogeneous distributions of 
precipitation and population and because depletion of the cloud by a rainout 


event in one part of the threat area reduces the threat elsewhere.) 


The roles of debris transport and dispersal are subtle, and one often 
eannot know in advance what condition is or is not conservative because complex 
tradecffs determine rainout effects. For example, higher wind speeds drive the 
debris cloud over more people, but the hazard to each is less because of the 
reduced cloud residence time overhead. Lower wind speeds, on the other hand, 
imply that the cloud will threaten fewer people but will present a greater 
hazard to those people threatened before it is dispersed. Some certain wind 
speed will pose a maximum danger, but what that speed is depends on other 
driving factors (yield, dispersal rate, and population distribution) and can be 
determined for a specifie scenario only by sensitivity and tradeoff studies. 
Similar comments apply to dispersal rates: the more rapidly a cloud grows, the 
more area it will cover, but the more rapidly it will be diluted, and vice 
versa. Dispersal rate itself also depends on the prevailing wind field pecause 
of shear effects. 

The meteorological factors that determine debris dispersal and transport 
are, respectively, small- to mesoscale atmospheric turbulence and meso- to 
large-scale wind fields, the latter representing a modeling area where 
appropriate climatic data can and should be introduced. 

C,. ‘Time and Space Distributions of Debris-Precipitation Interactions. These 
topics include most of the meteorological aspects of the rainout problem. we 
ean further separate the physical elements most important to rainout into two 
interrelated but distinctively different groups: (1) factors that determine 
the frequency, location, and extent of interception regimes between debris 
clouds and relevant precipitation cells, and (2) influences of local 
circulation features in the vicinity of debris-cell encounters (for example, 
entrainment or detrainment). Group 1 dominates because it directly determines 
the likelihood, locaticn, and overall coverage of rainout events. Group 2 
arfects the intensity of a rainout event and, to a lesser degree than group 1, 
the extent and location of ground deposition. 

Coincidence probabilities and intersection regimes for debris and 
scavenging-cell encounters depend on the mesoscale time and space structure of 
precipitation cells. This structure depends in turn on position within a 
synoptic~scale disturbance (Appendix A) and ultimately, therefore, on regional 
climate. Knowledge of the time-space trace of a debris cloud and of the 


stochastic extent and behavior of precipitation cells (from precipitation- 


morphology data and modeling studies) allows statistical assessment of the 
pertinent intersections. (The situation is not generally one where a small 
debris cloud is embedded in continuous or constant precipitation, but rather 
one where a chance encounter occurs between two widely varying but comparable 
volumes in space. ‘The complete superposition and mixing of these volumes is an 
improbable special case, so we must consider the extent of such intersections 
as well as their associated probabilities. ) 

Given an intersection betwecn debris and a precipitating system, it is 
generally also necessary to evaluate other aspects of the encounter. Air 
movements within and near the precipitation system will alter the local 
distribution of affected debris and should be accounted for in establishing the 
extent and intensity of resulting ground deposition. Changes in the 
precipitation system caused by the presence of the debris particulates, 
however, are not expected to be important. 

D. Seavenging Rates. Many microscopic physical processes exist by which 
precipitation scavenging could remove debris particles from air, and their 
relative magnitudes and roles are still poorly understood; nor are the 
pertinent properties of debris particles well established. These uncertainties 
translate directly into fundamental uncertainties in removal rates for the 
practical problem. 

Experimental work relies heavily on ground observations thrcugh broad 
hypotheses concerning how various removal mechanisms combine and how they vary 
with altitude, and how scavenging data for one particle species can be 
extrapolated to other, grossly different species. Unequivocal data are 
unlikely because airburst tests with nuclear devices are by treaty impossible. 

The existing literature indicates that the gap between theory and 
experiment also has not been bridged. The large variety of possible scavenging 
mechanisms and their dependence on medium- and small-scale variations in 
moisture, particle size, particle physical-chemical characteristics, 
temperature, airflow, and storm-cloud life cycle probably mean that theory is 
many years away from providing a definitive substitute for real airburst data. 
Indeed, theory is probably several years away from synthesizing the separate 
iachantsns into a form we can directly compare with present field observations. 

Detailed studies by other rainout investigators3 imply that debris 


scavenging will occur with efficiencies that are in the range of maximum 
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rainout danger, with an uncertainty factor of %10 in absolute scavenging 
efficiency. Tne rainout problem therefore cannot be dismissed solely on the 
basis of arguments that scavenging rates are too low for airburst debris to 
pose a problem, at least until a more complete and convincing synthesis of the 
whole scavenging-rate problem becomes available. 

we can bound rainout effects, however, through sensitivity studies made 
with our assessment system; and our results show that collateral-damage effects 
are relatively insensitive to overall scavenging efficiency over the 
“dangerous” range. The whole problem can then be put in proper perspective 
with respect to the other driving elements. 
EB. Population Distribution. The importance of local population density and 
distribution was pointed out earlier. whe need a detailed, realistic, 


population data base to estimate rainout collateral-damage effects rather than 


threat areas. 


EF. Postdeposition Factors. kven if we know the amount and distribution of 
debris deposited at the earth's surface, we still do not have a direct measure 
of real doses that resident populations would receive. We should account for 
several other phenomena, most of wnich would be protective in character. 

For example, runoff and washoff from hard surfaces would transport 
debris-laden precipitation to storm sewers and rivers, thereby reducing 
casualties in just those areas with the most concentrated populations, the 
cities. Similar comments apply to shielding by structures. Ground absorption 
would be effective in reducing dose rates in other areas, as would shielding by 
local terrain features. 

Another broad area of importance is protective response by threatened 
people. For example, the potentially.longer delay times inherent with rainout, 
as compared with those for prompt effects or dry fallout, mean that significant 
time (perhaps many hours) would often be available for protective measures 
before large doses were accumulated. # 

Proper assessment of relevant postdeposition factors is a major task that 
is beyond the scope of our present rainout study. Wwe can make such assessments 
later, however, because they are essentially independent of the other elements 
in our system; and, because the overall effect would be to reduce rainout 
collateral-damage estimates, we have thus far viewed the whole area of 


*This time factor is another reason why rainout effects are comparable with 
prompt effects and dry-fallout effects only on a qualified basis. 
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postdeposition effects as a built-in conservatism ("safety" factor) in our 
results and have used the same conventional, unshielded measure that other 
investigators use in determining effects, the “infinite whole-body dose" (IWBD) 


(Appendix B). 


Vi. ThE ACRA SYSTEM 

Tne analytical tools we have developed to address rainout problems are 
embodied in a computer system we call ACKA, for "Atmospheric=Contaminant 
Hainout Assessment." The basic criteria for ACRA are as follows. 
Modeling Objectives: . 

e Emphasize pragmatic goals and solutions. 

@ Develop a broad, flexible, overall modeling and computational 
framework. 

@ Accommodate requirements for in-depth treatment of many difficult, 
interacting, physical elements. 

@ Eliminate as many irrelevant elements as possible by scoping their 
influences to make the problem tractable. 

Ground Rules: 

e Carefully define collateral-damage measures. 

@ take conservative, practical, state-of-the-art evaluations of driving 
elements. 

@ Rely on detailed work of others where possible. 

@ Use adequate and appropriate data bases (recognizing especially that 
time and space resolution of data must match the scales of pertinent 
physical events). 

@ &mphasize actual effects assessments and the implications of realistic 
operational contexts. 


Operational Approach: 
@ Use scoping and bounding calculations to identify potential 


collateral-damage situations. 
@ Make finer-cut sensitivity studies to discriminate and evaluate 
systematically the roles of individual driving elements. 
e Use information feedback to upgrade model components systematically as 
necessary to improve the quality of assessment results. 
.@ Iterate until the "problem" goes away or the results become 


ag. 


definitive. a 
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The ACkKA system (Fig. 2) consists of eight computer codes, each composed 
of easily replaced modules. These eight codes are independent but communicate 
with each other (Sec. VII) to give an integrated overall system. The eight 
codes and thei:* functions are as follows. 

GEBESG: Prepares input-data libraries for SpPOOr. 

SPGUK: bebris transport and dispersal simulator. 

HESPOUOR: eprocesses SPUUK output for input to SCAWUEN,. 

TEMPEST: Storm and precipitation-cell simulator. 

SCAWUEN: Scavenging simulator. 

SHASTA: Tallies rainout effects from a single shot.*# 

SHASTAX: Tallies rainout effects from multiple shots. 


RAIN: Tallies rain amounts at fixed ground points. 


*The codes GEbéSU through ShASTA calculate a single input source configuration. 
Normally, we use a source representing the stabilized cloud from a single shot. 
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A. _GEBESQ. GEBESO uses input atmospheric-turbulence spectral data to prepare 
autocovariance functions and corresponding Bayi parameters for the debris- 
particle random-walk recursion relationships in SPOOR. These parameters have 
been described elsewhere.4 

GEBESO extends the original SPOOK code to permit an arbitary number of up 
to 100 correlated velocity terms in the recursion relationships.® It also 
provides an automatic capability for building libraries of Be sets (for 
different input spectra and various numbers of correlation terms) that can be 
selected by SPOOR. 

The input to GEBESO is an atmospheric turbulence-energy spectrum, and its 
output is a set of Bed values for SPOOR. 
Bb. SPOOK. This code is described in Ref. 4; only an outline will be given 
here. SPOGK represents an initial source cloud with a large number of 
individual particles. Each particle is "tracked" by simulating its wind-driven 
trajectory through a series of small time steps to produce a set of independent 
random-walk trajectories in three-dimensional time and space. ‘This tracking is 
done by statistically sampling simulated local wind and turbulence fields as 
each particle advances in time. These fields are based on generalized wind 
data for large-scale flow and on turbulent-energy spectra for micro- and 
mesoscale flow. To include effects from a broad spectrum of turbulent-energy 
scales, input for the large-scale field is derived from a climatological data 
base for Western Europe (Appendix C), and micro- and mesoscale wind 
fluctuations are introduced through a power spectral density. he thus simulate 
the role of turbulence in both meander and dispersal, naturally account for 
complex flow fields and time-dependent diffusion rates, and automatically 
simulate three-dimensional shear effects in the ensemble of particle 
trajectories. 

This model forms the core of the ACRA modeling approach. Not only is it 
inherently realistic, but, more important, it also makes tractable the 
mathematical problems of effecting the probability combinations (Secs. III and 


IV). The model generates independent random-walk trajectories for the 


*Only three terms were allowed in the original SPOOK code. 
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Simulated debris particles and implicit statistical properties of the ensemble 
of trajectories thus replace the explicit “cloud" properties of other 
formulations. because each particle traces out a well-defined, pointwise, 
history in time and space it has a straightforward relationship to other 
variable fields defined in the same phase space, such as precipitation-cell 
morophology, precipitation intensity, scavenging probability, and population. 
There is no need to resort to convenient but artificial debris-geometry 
representations as a basis for estimating interactions. This is in contrast 
with the more conventional deterministic modeling frameworks, where simulation 
of stochastic interactions between a "cloud" and the other problem variables 
becomes difficult and somewhat arbitrary. 

In simulating material transport in any domain that is subject to a series 
of interdependent random processes, the tracking, interaction, and tallying 
techniques must also account for the esSential stochastic variabilities without 
biasing the result by intermediate integrations, such as defining a "cloud 
size" before accounting for precipitation intercepts. It is intrinsically 
realistic to generate a particle trajectory along which one or more of these 
processes can occur and then impose the proper phase=-space-dependent process 
simulations and build up the "cloud" effects as an ensemble of these individual 
elements. 

The input to SPOOK consists of a set of Baya values from GEBESO, a set of 
small-scale wind-fluctuation variances, a large-scale wind data base, and the 
seenario data: yield, stabilized source-cloud dimensions, source-cloud 
altitude, initial debris distribution, shot time, and shot geographical 
location. Outputs from SPOOR consist of state variables (time, position, 
velocity) for all particle trajectories at all time steps, for input to 
RESPOOR, and a summary printout of the first and second moments of the state 
variables at each time step. 

C. AkSPOUR. HeESPOUR is a data-handling and processing code. Its function is 
to reorder and sort the SPOOK output data to permit more efficient processing 
by SCAWUEN. 

D. TEMPEST. TEMPEST is a dynamic, Monte Carlo, storm- and rain-cell 


morphology simulator. It sets up a fixed, user-specified storm geometry 
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composed of any number of separate "zones"* and then generates a characteristic 
set of moving, elliptical rain cells within each. These cells are 

transported with inputespecified velocities within each zone, and they start, 
grow, and die out according to input PDFs for initial cell size, cell aspect 
ratio, cell orientation, maximum cell size during its lifetime, and starting 
life-cycle status of the cell. Appendix D gives specifics of the TEMPEST 
model. 

1. Coordinate Systems. because different parts of our code system came 
from different sources, we have used three coordinate systems to avoid major 
recoding efforts, | 

The basic geographic (I,J) coordinate system and a corresponding 
rectilinear system (x,y) are shown in Fig. 3; the (I,J) coordinate system and 
the (1,J})-to-(x,y) transformations are described in Appendix C. 

The SPOOR, SCAWUEN, RAIN, SHASTA, and SHASTAX codes use (I,J) and (x,y) 
coordinates, ** TEMPEST, on the other hand, uses its own rectilinear 
coordinates (X,Y¥}, also shown in Fig. 3; all TEMPEST calculations are performed 
with respect to that system. The system origin is defined by the storm- 
geometry inputs. The transformation between (x,y) and (X,Y) coordinates is X = 
“-y, Y= xX. 


a i cn ne ee a 


*#we have incorporated une geometry package from the MCN® family of LASL Mor.te 
Carlo codes into TEMPEST. (This sophisticated peqneery system allows us to 
define any configuration of storm zones described by quadratic surfaces. 


*#Transformations from (LAT, LONG) geographic coordinates to (I J) are also 
included in some instances, such as for SPOOR scenario inputs. 


X ("NORTH") 
Fen J, x, OR Y ("EAST") 
I OR y 
Fig. 3. 
The (I,J), (x,y), and (X,¥) coordinate 
systems. 
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2. Storm Geometry Example. Figure 4 shows a storm-system geometry that 
characterizes the important features of the composite storm structure of 
rig. A-1. Wwe used a simple four-ellipse configuration to construct the storm 
zones in Fig. 4, although, for general purposes, the geometry package in ACRA 
offers much more flexibility. 

The synoptic storm center is as indicated near the intersection of zones 
1, 5, 4. zone 1 represents the ragged stratiform precipitation region, whereas 
zone 2 includes the area of light, continuous precipitation discussed in 
Appendix A. Zones 3 and 4 respectively depict the cold frontal band and the 
postfrontal shower region. 

tigure 5 shows the same storm with an overlaid set of representative 
rain-cell centers and typical cell velocities. Note the characteristic 
zone=to-zone variations in rain-cell density. The number of cells is 
determined by using the methods described in Appendix D and the concepts 
outlined in Appendix A. Zone 1, the ragged stratiform region close to the 
storm center, has the most dense precipitation and is where the heaviest 
"continuous" rains usually fall. Zone 4, the postfrontal region, typically has 
the least rainfall coverage. The coverage in this region (as in the other 
regions) varies greatly from storm to storm but is perhaps more dependent on 
local geographic and seasonal conditions than is coverage in the other regions. 
Cold air flowing over relatively warm, moist regions will often cause 
more extensive precipitation coverage in zone 4, both in area and density. 
Northern Germany, with winds fetching moisture off the North Sea, may have 


relatively high winter precipitation coverage under northwest flow (Zone 4). 


Fig. 4, Fig. 5, 
Tllustrative storm geometry. Representative rain-cell centers and 
zonal cell velocities. 
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A storm structure such as that shown in Fig. 5 can be used to demonstrate 
the role of climatic storm tracks in determining precipitation behavior at a 
particular site. For instance, Woronicz® shows primary and secondary storm 
.tracks for Western Europe for each month of the year. All months are 
characterized by two preferred tracks in the area of Uermany. A northern storm 
track exists from Scotland through southern Scandinavia, and a southern track 
runs across rrance and italy. The nortnern track would subject German stations 
to a preponderance of transections across Zones 3 and 4. The southern track 
might put German sites in zone 2, but more often they would be in the far 
regions of perhaps a more extensive Zone 4, where onshore northwest winds might 
bring weak, shallow showers. 

The cell movement reflects a circulation about the storm center 
representative of middle tropospheric winds. Appendix A suggests that this 
relationship between cell movement and mean airflow is simple and useful. 

3. Mesoscale Hain-Cell Structure. Our goal with the TEMPEST simulator is 
to provide a dynamic structure of rain-cells that debris particles must 
traverse, a structure whose statistical properties match what we know about 
such structures in nature and whose basic parameters we can adjust to allow 
normalization of the overall ACRA system to pertinent irtegral data, such as 
average rain-area fractions, rainfall frequency distributions, and average rain 
rates. 

A main problem, therefore, is to specify statistically the time-dependent 
geometric properties of this structure (cell position, size, number density, 
and shape with respect to the mean flow). A concomitant problem is to specify 
how the properties of the precipitation field depend on position in the 
large-scale synoptic storm pattern. These topics, discussed in Appendix A, 
form the bases of the TEMPEST model. 

TEMPEST provides a field of statistically generated precipitation cells in 
the same phase-space domain as that traversed by the SPOOK trajectories. These 
two fields can be combined to establish locations and frequencies for 
debris~-particle encounters with precipitation cells of appropriately varying 
size, density, and intensity. 

. 4a TEMPEST Inputs and Outputs. The input for TEMPEST consists of the 


following. 
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@ Storm-gecmetry specifications. 

@ kain-cell velocity for each storm zone. 

@ pbs for rain-cell aspect ratios, maximum céll size, and cell 
orientation with respect to cell velocity.. 

TEMPEST output consists of the state variables (position, size, 
orientation, ana life-cycle status) of all rain cells at all time steps for 
input to SCAwUEh, plus printouts of selected cell summaries for monitoring’ ° 
purposes. 
be.  SCAWUEL 

This code transforms all the TEMPEST cell data into a moving storm system 
in geographic (SFOUR) coordinates (Fig. C-2, for example) and then assimilates 
the KkKnSkbGUK particle data, identifies resultant interaction sites for debris 
particles in precipitation cells, calculates local within-cell scavenging 
probabilities, produces a corresponding set of scavenged-particle weights, 
transports these scavenged particles to the ground, tabulates their statistical 
properties, sets up a Lagrangian tally grid based cn these properties, and, 
finally, tallies the intensity of ground-deposited debris on this grid. These 
functions of Scawuen are more clearly understood in terms of the logic-flow 
diagram of trig. 6 and the specific SCAWUEN models described in Appendix &. 

lnput for SCAWUEN consists of the shobetine location on the TEMPEST time 
axis, the starting position and velocity of the TEMPEST storm in SCAWUEN 
coordinates, rain-intensity parameters Ye and Vy vertical-intensity tables 
a(4), a horizontal velocity Vos for transport of scavenged particles, PDr 
tables for the maximum cell-intensity parameters Hos and the tally-bin 
parameters k, NW, and rt’. The output of SCAWJEN is the dimensions, geographic 
locations, and deposited-debris intensities of a set of ground tally bins for 
input to ShAS’A, plus corresponding printouts of user-selected portions of the 
intermediate and final SCAWJEN results. 

Be. SHASTA, SHASTA tallies collateral-damage effects. It uses the output from 
SCawUEN and an input population data base to do three things: (1) identify and 
locate population points within the ground-deposition tally grids, (2) 
ealculate “infinite whole-body dose" (Iw5D) contributions for all time steps 
and each affected town, and (3) find the total IkbbD for each town, summed over 


all time steps. 
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Fig. 6. 
Functional flow diagram for SCAWUEN. 


The population data we use are taken from RAND's TANDEM data base for 
Central Kurope. ! These date are incorporated into a specially prepared 
population data base that we designed specifically to minimize computing time 
in SHASTA. 

ShASTA estimates IWEL (see Appendix b) as 

d 


jk IwbD at population point k from deposition during time step j 


4 ~0.24 ¥ 5 
2.02 x 10't, uyk ? 


where Th = deposited-detris intensity of the tally cell within which popu- 


lation point k is located at time step j, 


XY = total fission-equivalent yield of the source cloud, 


M = number of SPOOK particles in the calculation, 
.t; = time from fission. 


The total Llwbb at point k from the source cloud for all time steps is then 
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where N is the total number of time steps in the SPOOK simulation. 

Input to SHASTA is the ground-debris information from SCAWUELN and the 
ShASTA population data base. The ShASTA output is a set of data summaries for 
eacn affected population point; these summaries include geographic location, 
population, and b, for each town, sorted and listed to facilitate follow-on 
user analyses. A final tally lists the cumulative rainout-affected population 
totals within specific Dy contours. 

G. SHASTAX. This code combines the output from multiple ShASTA calculations 
to give total effects tallies for multishot simulations. SHASTAX makes linear 
superpositions of the individual D, levels from an arbitrary number of SHASTA 
runs to tally total multishot Iwbb for affected population centers. ‘the output 
format of SHASTAX is similar to that of SHASTA. 

Hh. KAlN. We designed our storm and rain models to be fundamentally literal 
and te include certain basic, adjustable parameters so we could use a broad 
spectrum of meteorological data (1) to verify the models and establish proper 
operating regimes for our system parameters and (2) to determine the real-world 
climate-based probabilities we can attach to our simulation and sensitivity 
studies (see Sec. VII1). 

the RKA1LN code addresses botn of these goals. The specific function of 
HAIN is to extract explicitly the time history of ground rainfall implicit in 
any TEMPEST-SCAWUEN simulation. Ground-rain data such as frequency 
distributions and time-averaged rates are implicit statistical properties of an 
ALKA simulation just as they are in a real-world situation, and these 
properties must be derived or "measured" by Simulating the real-world 
experiments that are used to generate the meteorological data.* 

ln the KAIN code, we tally total rain at a fixed ground point ("rain 
gauge") during each time interval At in a TEMPEST-SCAWUEN calculation. RAIN is 
essentially the SCAWUEN code with modified input and tallying modules; it 
substitutes the input ground point for the normal KESPGOR output, uses a 
standard TEMPEST output, and, fer each At, time-integrates (over At) the rain 


*#The same statistical and the same time- and space-resolution considerations 
also prevail in the ACKA "measurements" as in the real world. 
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intensity from all cells passing over the point during At. Appendix F 
describes details of these integrations. 

The KAIN output is a listing of total rain at the ground point during each 
At. We can bin and average these data as desired to provide frequency and 
average-rain-rate information with any desired time resolution between the 
limits At-and T (total simulation time).# 

VII. ACKA OPERATIONAL FEATURES 

The eight independent ACKA codes are designed so that communication 
between appropriate codes is achieved through conjoined input-output 
designations. Figure 7 shows the ACKA functional flow diagram. In the general 
case, all codes may be exercised in one action to obtain end results. ACKA 
does, however, give the user the flexibility to “save" outputs from any code 
and terminate the execution stream at that point so that he can review 
computational results before executing the next code in tne sequence. 

A typical course of action for a user studying rainout effects with ACKA 
might oe first to use the TEMPEST code to generate a number of storms that have 
precipitation eells of various desired cheracteristics. After reviewing the 
results, the user might make several runs with the RAIN code, varying 
rain-intensity inputs and using any of the sturms previously generated. Once 
the user feels he has properly modeled the storms, he is ready to continue his 
analysis. 

After identifying the appropriate atmospheric-turbulence spectral data to 
be used as input, the user would make a G&btsU run to generate the By yi 
Parameters needed as inputs to the SPUUH code (debris-clcud simulator). At 
this point, he might use SPOUh to generate a set cf debris clouds and then make 
a XESPCOR run with the SPUOUK outputs to build input files for use with SCAWUEN. 

Using any of the debris clouds and any of the storm systems, the user can 
then execute the SCAWUEN code (scavenging simulator) and the ShASTA code 
(tallies rainout effects from a single shot). A final user option is to 
execute the SHASTAX code to tally rainout effects from multiple debris clouds 
by combining, in any order, the ShASTA results obtained for the single shots. 

we nave optimized all eignt programs to maximize executional efficiency 


and reduce overall execution time (see Table 1). 


#\ie typically use At = 0.1 h and T = 5-10 h. 
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Fig. 7. 
ACRA functional flow. 
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TABLE I 


EXPECTED RANGE OF EXECUTION TIMES FOR CODES 
IN THE ACRA SYSTEM 


(CDC 7600) 
Range of Execution Time 
= min 

Code Name Minimum Maximum 

GEBESO » ae 5.0 

SPOGH 0.9 4,0 
‘ RESPOOR 2.0 8.0 

TEMPEST 1.0 3.0 

RAIN 0.4 1.5 

SHAWUEN 0.5 3.0 

SHASTA 0.6 2.0 

SHASTAX 0.5 3.0 


VIII. WAINOQUT CALCULATIONS~-SCOPING AND SENSITIVITY STUDIES 

Because a bewildering number of physical variables influence the 
real-world collateral-damage effects from rainout, an exhaustive examination to 
determine the differential roles of all these variables for rainout assessments 
would be impossible. Fortunately, however, rainout collateral damage is an 
inherently integral phenomenon and therefore is relatively insensitive to 
details as to how most of the pertinent physical variables actually vary with 
time and space, at least over regimes of practical military interest. 
Moreover, even if we knew the exact sensitivities of rainout effects to these 
variables, they would generally be masked by uncertainties in the statistical 
fluctuations introduced by the fundamental meteorological variabilities 
(See. IIT). 

we do not mean to imply that one can simply dismiss such factors from 
eonsideration, but rather that most of them will, in fact, be of secondary 
importance in rairiout assessment, provided we properly understand and scope 
their roles. Wwe believe this scoping to be an important task, as evidenced by 
the substantial effort we have expanded in scoping studies and by the 


complexity of the assessment system we evolved. 
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A major part of our effort has therefore been directed to progressively 
eliminating as many irrelevant variables as possible and to identifying the 
driving elements and perametrically evaluating their roles in raincut-effects 
assessment. Certainly, we have not been totally successful in achieving these 
somewhat idealistic goals, and several areas may require more work. 
Nevertheless, we have made major progress toward these objectives in developing 
and implementing the ACRA modeling concepts, and both the overall character and 
the details of this comprehensive calculational system reflect a host of 
insights we have gained during our studies. We have also used the ACRA system 
in its present form to simulate detailed rainout collateral-damage effects for 
over 500 representative nuclear shots to examine how these effects depend on a 
variety of basic physical and modeling variables. we will discuss the results 


of these simulations in the following sections. 


A. Single-shot Studies. 
1. A Typical ACKA Calculation. frigure 8 shows the "threatened" 


population points for a 1-kt fission shot in west Germany (shot No. 2 of the 
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Population points within 40 km over which debris passes; 
single 1-kt shot. 


referenced scenario discussed in Sec. VIII-b-1). The figure indicates the 
initial stabilized source cloud, the path followed by the centroid o”% the 
debris cloud, and those population points within 40 km over which some parts of 
the cloud passes.* the mean wind speed for this case is 35 km/h, in the 
direction along the centroid path. Table III (p. 40) gives the source-cloud 
dimensions. fopulation levels for the 54 indicated points range from 31 to 
4400; the total population for all 54 points is 44% 471, (No debris cloud is 
shown in the figure because our debris clouds consist of three-dimensional 
time- and space-varying ensembles of individual points.) This tally of 
"shadowed" population points is effected by bypassing the rain-cell and 
scavenging simulation parts of TEMPEST and SCAWUEN, through simple adjustments 
of input parameters. 

Figure 9 shows a population-point plot for the same 1-kt shot used in 
Fig. 8, but now with a typical TkMpeST-SCAWUEN rainout calculation included to 
Simulate a medium-rain situation for west Germany. the designations light, 
medium, and heavy rain are used to indicate different rain area fractions (and, 
hence, different average rain rates). tigure 9 shows only those population 
points that receive rainout deposition, again for only the first 40 km.** The 
12 points in Fig. 9 range in population from 156 to 1400; the 12-point total is 
6042. Figure 10 plots the number of people affected in the calculation as a 
function of infinite whole-body dose. 4 single, "unlucky," high-dose event is 
indicated at D = 520, N = 700; this event occurred at the point indicated by 
the arrow in Fig. 9. 

Another useful type of plot (rig. 11) shows dose versus distance from the 
shot point for the population points affected by rainout.*** This figure 
illustrates four characteristics of any rainout situation: (1) affected points 
have a bounding "envelope" that is determined primarily by yield and mean wind 
field; (2) the envelope is very nonlinear at distances within a few tenths of 


an hour of cloud-travel time from the shot and has a characteristic range of 


*For convenience, we limit the plot in rig. & to 40 km. ‘he calculation, of 
course, extends well beyond that limit. 


**¥Four points in Fig. y are different from those of Fig. &, because the rain- 
rate-dependent, scavenging-site-to-ground transport included in the rainout 
calculation causes a small change in the distribution of ground-deposited 
debris. The population points affected in the two cases are thus not guite 
identical. : 


##*This is the same calculation as for figs. 9 and 10. 
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Population points within 40 km that receive rainout deposi- 
tion; single 1-kt shot in medium rain field. 
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Dose-vs-population distribution for illustrative calcu- 
lation: single 1l-kt shot in medium rain field. 
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Dose vs distance for illustrative cal- 
culation; single 1-kt shot in medium 
rain field. 


maximum IwWBD threat that is determined primarily by yield and initial cloud 
size; (3) there is a maximum range beyond which dose falls off very rapidly 
because of depletion and cleud dispersion; and (4) most populetion points 
receiving rainout doses are weil within the bounding enveélcpe. 

2. Statistical and Numerical-iieshing Studies. Several of our early 
code-development and single-shot-sensitivity studies concentrated on tradeoffs 
between computer running time and numerical-meshing and statistical~sampling 
parameters. This work established meshing schemes and sample sizes required to 
achieve adequate numerical accuracy within practical computing times. We 
examined the effects of the following: 

@ Number of 8 terms in the SPOCR recursion relationships. 

Number of SPOOR source particles. 
Tally-bin radial-meshing models. 


® 
e 
@. Tally-bin mesh sizes. 
@ Time mesh size. 
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Fig. 12. 
Effect on rainout results of varying source sample size. 


We found that rainout results and computer running times were insensitive 
to the number of 8 terms in SPOOR, for numbers between 3 and 40. We used 10 
Bts in most of our subsequent work. 

Computing times in the ACHA system vary in a roughly linear fashion with 
the number of source particles M, as we would expect. Figure 12* shows results 
for a set of point-source calculations designed to study the effect on rainout 
results of varying M. Wwe found a few thousand particles to be generally 
adequate; most of our simulations use Ms 5000. 

Developing an adequate tally-bin scheme is difficult because extreme 
debris-concentration gradients in both time and space must be treated to 
achieve adequate resolution in the results, particularly near the origin. A 
Lagrangian tally grid is a necessity; even then, one is still trying to resolve 
a very complex gound-deposition field whose spatial distribution is essentially 
unknown. We found the second moment of ground-deposited debris about the 
centroid to be useful in characterizing these spatial distributions, and, after - 


experimenting with several semiempirical meshing models, we settled on the one 


*uur sensitivity results are often presented only in terms of relative dose 
effects in arbitrary units, when only relative results are available cr per- 
tinent. Absclute rainout effects are presented in terms of tissue doses (rem). 
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described in Appendix E. With this scheme, rainout results appear to be 
relatively insensitive to tally-bin size. This insensitivity is shown in 

Figs. 13, 14, and 15 for 1-, Z-, and 10-kt shots, respectively. An inherent 
pitfall in using small radial-mesh sizes is shown in the 40x(40+20) curve of 
Fig. 15. Near the shot point, statistical fluctuations in the rapidly varying 
Spatial gradients of ground-deposited debris can give abnormally small oe with 
correspondingly small bin sizes and aberrantly large dose estimates (see 

Eqs. B&-13, B-14, and E-15). This phenomenon is strictly a statistical anomaly 
and is inherent in the tally-bin meshing scheme. These errors are significant 
only within the first few tenths of an hour from the origin, but the effect 
normally causes us to restrict the number of radial bins to no more than 145-30. 
fmpirically, this range gives a good compromise between overmeshing the cloud 
near the origin and undermeshing it elsewhere. 

Another source of uncertainty in our calculations, again within 0.3 h of 
cloud-travel time from the shot point, is in the time-mesh size, At. We 
typically use At = 0.1 h to avoid the larger computing times required for 
smaller At's. As the calculation cannot give good resolution for the first two 
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Effect on rainout results of varying the tally mesh; single 
l-kt shot. 
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Fig. 14. 
Effect on rainout results of varying the tally mesh; single 
2~kt shot. 
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Fig. 15. 
Effect on rainout results of varying the tally mesh; single 
10-kt shot. 


or three At's, there is an uncertainty in aetailed dose calculations over this 
range, the same range discussed in the previous paragraph. A solution to this 
problem, as well as that in the previous paragraph, would be to use a finer At 
near the origin; that is, use a time-varying time mesh. This option is not 
presently available in ACRA, but it would be a useful extension of the systen. 

we. have made empirical studies of the two uncertainties mentioned in the 
previous paragraphs by using different At's and a variety of tally-bin meshes. 
These studies have led us to consider our dose results for the first few tenths 
of an hour from the origin to be uncertain by a factor of up to 2 or 3 
(generaiiy on the high side.) 
3. Initial Sensitivity Studies. Several of our early single-shot series 
‘ also examined physical factors that we expected to be of secondary importance 
for determining rainout effects. These factors included the following: 

@ Spatial distribution of activity in the source cloud. 

Input turbulence spectrum for SPOOR. 
Velocity of the TEMPEST storm. 
Velocity of TEMPEST rain cells. 


Air-to-ground transport of scavenged debris. 

we used various combinations of uniform and Gaussian activity 
distributions, in both the axiul and radial dimensions of representative source 
clouds, to examine the effects of source internal distribution on ground doses. 
These effects were small, because the initial activity distributions 
disappeared so quickly as a result of turbulent diffusion that a debris cloud's 
"memory" of its initial internal configuration lasted only a small fraction of 
an hour. Initial cloud size is, of course, more Significant. 

The role of input turbulence spectrum is less clear-cut. Higher 
turbulence spreads the debris cloud more quickly over more area than does lower 
turbulence, but the cloud density is also reduced more quickly. The net effect 
of these competing factors depends on the local population distribution; 
nowever, the differences resulting from different input turbulence spectra 
usually are not large and are masked by fluctuations in the other variables, 
notably population distribution. Atmospheric shear from vertical var.sation in 
the horizontal wind fields is more important and is automatically included in 


our calculations (see Ref. 4). 
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Similar comments apply to the importance of TEMPEST storm and rain-cell 
velocities.. Rainout effects appear to be relatively insensitive to reasonable 
changes in these variables. This point needs more study but thus far we have 
found no credible combinations of these variables that produce large 
sensitivities in the final rainout dose estimates. 

Figure 16 illustrates how including the air-to-ground transport of 
debris-laden precipitation affects rainout doses. Again, this factor seems of 
secondary importance, but we include it in our rainout calculations (Appendix 
B). 

4, Hain Parameters and Scavenging-Efficiency Studies. Our early studies 
led us to conclude that rain and scavenging parameters were of more basic 
importance than those discussed above for determining and characterizing 
rainout-effects threats. OUnce again, single-shot parameter studies with ACRA 
were useful for delineating and scoping the roles of these elements. The 
elements studied were as follows: 

@ Probability density function for Ky = maximum rain intensity at 

the center of a rain cell. 

®@ Yy, Yy = within-cell rain-intensity distribution parameters. 

@ scavenging coefficient. 

@ Yield. 

figure 17 gives the standard Ky distribution we used for most of our 
ealeculations. Figure 18 shows the effect of varying K, for a 0.4-kt shot in 
medium rain, and tig. 19 shows results for a 0.4-kt shot in neavy rain. 
figures 20 and 21 give analogous results for a 1-kt shot. These figures show 
several examples of relatively large steps produced by single population 
points, and such major steps must be considered specifically in making 
quantitative comparisons. The most notable example is evident at 1-3 rem, 
where a single point has a population of 200 000; Table I1 gives dose data for 
this point from the calculations in Figs. 18-21. When we examine the results 
in Figs. 15-21 in detail, we find that threefold variations in R, generally 
produce much less than threefold variations in IwbD. We conclude that Ro is 
not a dominating parameter in determining rainout effects. 

Figure 22 shows how varying the rain-intensity distribution parameter 
affects rainout doses for a 1=kt shot in a heavy-rain field, and Fig. 23 shows 


analogous results for a medium-rain zone. The parameter Y and the rain-cell 


33 


de Ty 


{-WITH TRANSPORT 
foe I 


WITHOUT 
TRANSPORT 


o 


Y (arb units) 


io 10) ae | 
10) 10° io” io" 10° IC 
POPULATION WITH IWBD 2 Y 
Fig. 16. 


Effect on rainout results of air-to-ground transport of debris- 
laden rain; single 0.4-kt shot. 
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Effect on rainout results of varying Ri Single 0.4-kt shot 
medium rain field. 
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Fig. 19. 
Effect on rainout results of varying Ri single 0.4-kt shot in 


heavy rain field. 
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Effect on rainout results of ~:rying Ri single 1-kt shot in 
medium rain field. 


STANDARD.R, 


o| STANDARD Ry xlS= 7 
10 


STANDARD R,+ 2 | 


POPULATION WITH IWBD 2 D 


Fig. 21. 
Effect on rainout results of varying Roi single 1-kt shot in 
heavy rain field. 
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Fig. 22, 
Effect on rainout results of varying y; single 1-kt shot in 
heavy rain field. 
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Effect on rainout results of varying y; single 1-kt shot in 
medium rain field. 
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TABLE II 
IwBD FGR A SINGLE LARGE POINT (POPULAT1UN 200 000) 


IwkbD (rem) 


O.4—-kt Shot 1-kt Shot 
Standard Standard 
Re Ro + 2 Ry x 1.5 Ro fo #2 Ry x 1.5 
Medium Rain <1" <1 <1 1.3 <1 <1 
heavy fain <1 2.5 1.6 2.8 2.3 1.4 


5 


aWe arbitrarily eliminate STi points from the ACRA printouts whose total 


iwbD is <1 rem. 


2 


7 m 
because of their influence on the small-scale, rain-rate frequency 


size-distribution parameters o* and m are of fundamental importance, primarily 


distributions implicit in our model. hainout effects are relatively 


insensitive to modest changes in ae and m; effects vary with y as shown in 


rigs. 22 and 23. Small-scale (%0.1-n) frequency distributions, on the other 
nand, are sensitive to all three parameters. These parameters can thus be 
fixed to an adequate level of uncertainty for rainout calculations by comparing 
rain frequency distributions from the model with rain data, as outlined in 
Appendixes A and G. We concluded from these studies that reasonable values for 
the rain variables were 1/Y = 4, of = 2.5, and m = 0.47. Rainout-effects 
estimates are insensitive to modest variations of the parameters about these 
nominal values (see Sec. VIII-b-2). 

The scavenging coefficient A = khO+8 is discussed in Appendixes E and h. 
we examined the sensiti ‘ of rainout doses to scavenging efficiency for a 
1-kt shot in a moderate-rain field by varying the coefficient K over a range 
from 0.05 (0.5% probability of being scavenged in 0.1 h by a t=mm/h rain) to 
31.5 (95% probability of being scavenged in 0.1 h by a t=mm/h rain). Figure 24 
shows the results. kainout effects are insensitive to uncertainties in A’ over 
the most threatening range (1 < K < 30) because competing factors tend to 
cancel one another in determining overall population-versus-dose tallies. A 
larger A means that debris concentrations deposited from a cloud before it has 


been depleted will be larger than when ) is smaller; the concomitant faster 
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Fig. . 

Effect on rainout results Sao A; single l-kt shot in 

medium rain field. 
depletion, however, means that the threat lifetime and overall threat area are 
Smaller. Also, the probability that debris will be deposited on a population 
center is smaller when a few intense depositions secur (high A) than when a 
larger number of less intense events occur over a wider area (lower A). 
bventually, of course, as X is decreased, rainout doses begin to fall off 
because of the reduced deposition rates. The net result of all of these 
factors is a variation of rainout effects versus A that tends to have a broad 
maximum, with modest falloff as scavenging rates become so large that 
deposition is essentially “instantaneous” and with more rapid falloff as 
scavenging rates become small. (we aiso found this same behavior in an early 
study® using a predecessor of ACKA; in that study, all debris clouds were 
deposited uniformly with a predetermined, constant A.) 

Our rainout calculations to date have concentrated on small fission 
weapons (4 kt and less), but we have also made both single-shot and multishot 
scoping studies as a function of yield (fission devices only). Table III lists 
the source-cloud sizes and heights that we used. These dimensions are also 


shown in Fig. 25, along with a plot of the rain-intensity vertical weighting 
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TABLE III 


ASSUMED SGURCE-CLOUD HEIGHTS AND DIMENSIONS 


Yield Cloud conte: heats pu iecneds 
(kt) —__(km) (km) _ km 
0.4 2.4 0.6 1,4 
1.0 2.9 0.9 1.8 
2.0 3.5 1.z 2.2 
5.0 4.6 1.8 3.0 
10.0 6.2 2.4 3.4 
20.0 ' 9.1 -2 4.3 


INITIAL 
DEBRIS-CLOUDS \ > 


Skt 


ALTITUDE (km) 


0 Ol 02 03 04 05 06 O7 08 09 10 
a 


Fig. 25. 
* Initial debris-cloud sizes and heights and rain-intensity 
vertical weighting function a(Z). 


function a(z).* Figure 25 also shows the range of a within whicn each source 
cloud lies and within which most of the debris remains as the clouds are 
transported and dispersed. Rainout effects as a function of yield are again 
determined by tradeoffs among several competing factors, particularly higner 
debris concentration and larger threat areas for higher yields versus reduced 
rain rates and lower scavenging probabilities at the higher altitudes 
associated with higher-yield debris. Hain intensities above «7-6 km, for 
example, are expected to be too small to produce substantial debris scavenging, 
and we limited our calculations to 10 kt for this reason (see the indicated a 
for 20 kt in Fig. 25, for example). Single-shot calculations versus yield are 
not quantitatively definitive, because statistical fluctuations inherent in the 
different wind fields and rain cells encountered by the different debris 
clouds, plus the different local population distributions affected, obscure 
detailed comparisons; we thus show no quantitative single-shot comparisons. 
Two points are clear, however, from these studies: rainout effects tend to 
fall off for yields below ¥1 kt, and they tend to level off at 15-10 kt. We 
will explore these points more fully in connection with our multishot results. 
b. Initial Multishot Studies. A realistic, statistical simulation approach to 
rainout assessment is subject to statistical fluctuations as in nature, and 
these fluctuations lead to difficulties in interpreting the results of 
parameter studies, particularly for single-shot studies. Ideally, we would 
avoid these problems by simply running a statistically large number of shots as 
a unit, to sample effectively the variable fields (rain-cell and population 
fields, for example), and then make comparisons among the combined multishot 
results. Shot-to-shot fluctuations would then average out, and overall 
variations would become small in comparison with the effects being studied. 
Large numbers of shots require large computing times, however; so we have 
restricted most of our multishot calculations to 10-snot salvos as a useful 
compromise between increasing the computing time and decreasing the statistical 
uncertainties in the salvo results. 

1. An Illustrative Muitishot Calculation. Figure 26 shows the assumed 


initial source-cloud positions for our standard ten-shot scenario. These 


#See Appendixes & and I. 
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Fig. 26. 


Standard ten-shot scenario; cloud-centroid paths 
for 1-kt shots in typical wind field; rainout- 
affected population points within 40 km in medium 


rain field. 


locations were actually taken trom a simulated tactical nuclear engagement;* 


the only significance we attach to these locations, however, is that shot 


density and shot separation distance are probably typical for low-yield 


defensive nuclear engagements. 


Figure 26 also shows cloud-centroid paths for a 


salvo of 10 tekt shots in a typical wind field (mean wind speed of 35 km/h); 


population points within 40 km of the salvo that receive rainout deposition 


from a simulation in a light-to-moderate rain field (71 points with total 


population of “40 060); population points (circles) with IwBD > 100 rem (26 


points with total population of 117 000); and points (solid circles) with IWED 


> 450 rem (8 points with total population of 7 000), 


Figure 27 is the 


corresponding plot of population versus total IWBD. 
2. isuitishot Rain-Parameter and Seavenging Efficiency Studies. whe 


repeated several of our earlier single-shot sensitivity studies with multishot 


salvos to determine whether statistical uncertainties in the single~shot 
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*Shot 2 is located at LAT 50.7394, LUNG 9.6752. 
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Dose-vs-population distribution for illustrative calculation; 
ten l-kt shots in medium rain field. 


results were miSleading us in our earlier conclusions. Figures 2&8 and 29 show 
the effect of varying the p.obability density function for the maximum 
rain-intensity parameter Ky for 10-snot salvos of 0.4-kt and 1-kt devices, 
respectively; Figs. 30 and 41 snow the effect of varying tne within-cell 
rain-intensity parameter y for a series of 4-shot salvos; Fig. 32 compares 
results from 10-snot salvos for different rain-cell size-distribution 
parameters of and m; and tig. 33 shows the effect of varying the scavenging- 
coefficient multiplier kK. 

we found the same behavior in these multishot studies as we found in the 
earlier single-shot studies, and our previous conclusions thus remain 


unchanged. 


IX. MULTISHOT PRODUCTION KUNS. 

In the preceding sections we determired statistical and meshing parameters 
for ACKA, discussed the insensitivity of rainout effects to internal 
source-cloud details and to certain synoptic-scale dynamic storm features, and 


established reasonable, conservative, operating regimes for rain-morphology and 
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Effect on rainout results of varying Ri ten 0.4-kt shots. 
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Fig. 29. 
Effect on rainout results of varying Roi ten 1-kt shots. 
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Fig. 30. 
Effect on rainout results of varying y; four 1-kt shots in 
medium rain field. 
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Effect on rainout results of varying y; four 1-kt shots in 
heavy rain field. 
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Fig. 33. 
Effect on. rainout results of varying A; ten 1-kt shots. 


scavenging-efficiency parameters. Most important, we showed that rainout 
effects are relatively insensitive to variations in rain-structure and 
scavenging-rate parameters over the most threatening range and that large 
sensitivites appear only when the parameters are such that rainout effects 
begin to fall off rapidly. Wwe also demonstrated the fundamentally 
heterogeneous and probabilistic character of rainout collateral-damage effects, 
a natural consequence of the character of the driving rain-cell and population 
fields that determine where rainout depositions will occur. 

These studies and the ACRA system (exclusive of the RAIN code) deal with 
debris transport and dispersal, precipitation morphology, debris~cloud and 
rain-cell interaction regimes, scavenging and ground-deposition rates, and 
local population distributions. In the following sections we assess the 
remaining requirement for rainout, that there must be a precipitating storm 
system in the area. Wwe also examine the corollery question of the effect of 
being in different parts of a typical storm system. The goal here is to attach 
overall occurrence probabilities to our various rainout simulations through the 
use of regional climatic data. 

Gur approach is to calculate rainout effects as a function of input rain 
area fraction (HA) and then, using the RAIN code to estimate effective average 
rain rates versus RAF for our simulated storms, to interpret average rain rates 
in terms of the local climatic data to infer . robabilities of occurrence for 
given rainout situations. 

A. kainout Effects versus Rain Area Fraction. 

we again use our "standard" 10-shot scenario to estimate rainout effects 
for 0.4-kt and l-kt salvos. (higher yields will be examined in the next 
section.) These results are presented as a function of input KAF (to ITKMPEST). 
Actual rain area fraction is, of course, a time-dependent, statistically- 
varying quantity throughout any TLMre.Sl-SCAkUEN simulation. Table IV gives 
the input and actual rain area fractions for the present calculations. 

Figures 34 and 35 show rainout effects for 0.4-kt salvos when input RAfs 
are varied between 0.1 and 0.4 and between 0.015 and 0.2, respectively. figure 
34 exhibits a lack of sensitivity to RAF for values above 0.2; the 0.2-0.4 
RAPS correspond to a more than threefold variation in average rain rate, but 
the rainout-effects curves are all very similar (except for about a threefold 


variation in IwbD for a single population point 0.2 h downwind from the shot. 
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TABLE iV 


RAF lnrUT VALLES, ACTUAL RANGES, AND COHKESPONDING APPROXIMATE 
BPEECTIVE VALLES (RAE) FOR "PRODUCTION" RAINGUT RUNS 


Input Range of Actual HAFs Approximate 
HAF” for Production kuns Average (RAF) 
(0.910 0.009 to 0.012 0.011 
0.015 0.075 to 0.021 0.014 
0.03 0.029 to 0.031 0.030 
0.05 0.05 to 0.053 0.051 
0.10 0.12 to 0.13 0.12 
0.26 0.22 to 0.23 0.23 
0.30 0.354 to 0.35 0.35 
0.40 0.45 to 0.47 0.45 


line). sHKainout effects are somewhat lower when RAF = 0.05-0.1 and are falling 
off rapidly when KAF = 0.015. 

Figures 36 and 37 show analogous results for 1-kt salvos. The variation 
of these results with khAF is similar to that for 0.4 kt, although the rate of 
falloff at lower hkAF values is slower. 

These same trends can be seen more clearly in Taile V, which lists, for 
the calculations depicted in Figs. 34-37, the number of population points 
affectea at various iWbD levels. (The number of affected points is not a 
direct measure of collateral damage but is still a useful measure for 
identifying such comparative trends because it is not subject to the 
population-weighting fluctuations that are present in the final effects 
estimates. ) Tne data in Table V also suggest a broad maximum in rainout 
effects versus RAF at HAP = V0.1-0.2, for 1WKD of %100 rem and larger. This 
maximum is real, but it is masked in rigs. 34-37 by the point-to-point 
variations in population weignting factors that AURA attaches to the points of 
Table V to get the effects data in Figs. 34-37. 

The dependence of rainout effects on KaF seen in Figs. 34-37 and Table V 
is directly related to now debris-clcud sizes compare with distances between 


rain cells. Cloud size is essentially independent of KAr but does depend on 
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Effect on rainout results of varying input rain area fraction 
RAF; ten 0.4-kt shots. 
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Effect on rainout results of varying input rain area fraction 
RAF; ten 0.4~kt shots. 
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Effect on rainout results of varying input rain area fraction 
RAF; ten 1-kt shots. 
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Effect on rainout results of varying input rain area fraction 
RAF; ten 1-kt shots. 
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TABLE V 


NUMBLK OF POPULATION POINTS WHERE IWBD EQUALS OR EXCEEDS Do VALULS 
FOn LISTED RAFS, FOK 10~SHUT SALVOS OF FIGS. 34-37 


(0.4=-kt ShOTS) 


Do (rem) 

RAL 10 50 100 200 450 1000 2000 
0.015 18 2 4 0 ce) 0 0 
0.05 4G 26 19 10 6 2 0 
0.1 49 21 16 14 8 4 1 
0x2 54 13 1 9 5 4 2 
0.3 59 23 16 10 6 3 1 
0.4 78 27 14 6 3 1 1 

(1-kt SHUTS) 
De (rem) 
HAR 10 30 100 200 450 1000 2000 
0.01 25 9 4 0 0 0 0 
0.015 37 2 1 1 1 0 0 
0.03 23 le & 4 3 1 0 
0.05 9g 33 27 19 As) 5 1 
0.1 &5 34 25 20 14 5 ‘2 
0.2 150 42 22 12 6 2 0 
0.3 147 36 19 10 6 zZ 2 
0.4 174 4y 29 16 5 2 1 


yield, especially at early times, and increases rapidly with time; cell 
separation distance changes inversely with hAF and is essentially independent 
of time. Two characteristically different relative-size regimes thus result, 
with distinctively different dependences of rainout deposition on RAF. 

For RAF 20.1 and yields R1 kt, cell separation distances tend, on the 
average, to be significantly less tnan debris-cloud dimensions; and a deb:‘is 
cloud continually intercepts rain cells at a relatively high rate, so it "sees" 
a set of cells at all times that fairly well represents the overall rain-cell 


field. (In other words, the cloud is big enough to encompass a statistically 
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meaningful sample of cells from the large-scale field.) In this event, the 
properties of the ensemble of cells acting on the debris cloud will be 
Statistically "well behaved" (modest fluctuations) as functions of time and 
will vary with RAF in roughly the same way as the large-scale field. because 
deposition depends on cell intensity k only through the scavenging coefficient 
As KRO.8 | the KAF-dependent variations in the R of the intercepted cells are 
exactly equivalent to corresponding variations in 4; and the resultant 
dependence of deposition on i (aid, hence, on RAF) is weak (Secs. VILI-A-4 and 
V111-3-2), 

For RAF £0.05 and for lower yields, an early period exists during which 
rain-cell separation distances are larger than debris-cloud sizes and substan~ 
tial cloud dispersal takes place on the average before any rain cells are 
intercepted. Deposition thus begins at later times for lower RAFs because the 
cell separation distances are greater and because more time is required for the 
cloud to grow to a comparable size. Also, for smaller KAt's, the initial 
deposition rate is statistically more variable; and depletion, when it begins, 
is at a lower rate because of the lower average rain rate and the larger 
average cloud size at the onset of deposition. The net result with lower RAF 
is a rapid falloff of early-time deposition rates and an extension of the 
threat lifetime and, hence, of the distance and number of people threatened at 
lower LwWbb. 

Figures 34-37 and plots of LwWBb versus distance show these characteristic 
features for those population points affected by rainout (analogous to 
Fig. 11). Figures 38-40* show examples of such plots for the 0.4-kt salvos; 
note the differences in the maximum "envelopes." figure 41 shows analogous 


envelopes for the 1-kt salvos. 


BE. Fr" ut Effects versus Yield. 

‘Figures 42 and 45 show rainout effecte versus shot yield for "standard" 
10-shot salvos in a rain field with RAF = 0.10, and Figs. 44 and 45 give 
‘analogous results for a rain field with RAF = 0.015. Table Vl gives corre- 
sponding numbers of rainout-affected population points versus yield and IwBD. 


*In these figures, we plot IwbD versus distance from shot 2 for all points with 
IWbD > 10 ren. 
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Maximum dose-vs-distance envelopes for 


Dose vs distance for ten 0.4-kt shots 
ten 1-kt shots with various RAF. 


with RAF = 0.01. 


53 


10 10° 
FOPULATION WITH IWBD 2 D 


Fig. 42. 
Effect on rainout results of varying shot yield; ten-shot sal- 
vos with RAF = 0.1). 
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Fig. 43. 
Effect on rainout results of varying shot yield; ten-shot sal- 


vos with RAF = 0.10. 
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Fig. 44, 
Effect on rainout results of varying shot yield; ten-shot salvos 
with RAF = 0.015. 
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Effect on rainout results of varying shot yield; ten-shot salvos 
with RAF = 0.015. 
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TAbLE VI 


NUbEeR OF POPULATION POINTS WHERE IWBD EQUALS OK EXCERDS INDICATED 
D, VALUES FOR LISIED YIELDS (FOR 10-ShOT SALVOS OF FIGS. 42-45). 


(RAF = 0.10) 


YIELG D, (rem) 
{kt) 10 20_ 100 200 450 1000 2000 
0.4 4g 21 18 14 3 4 1 
1 &5 34 25 20 14 5 2 
177 47 36 22 15 8 1 
536 78 43 2g 17 9 5 
10 302 104 4g 31 14 7 3 
(RAF = 0.015) 
YILLD Dy (rem) 
(kt) 410_ 20_ 100 200 450 1000 2600 
0.4 18 2 1 0 0 0 0 
1 37 2 1 1 1 0 0 
2 107 14 7 4 1 1 0 
Y3 14 9 5 2 2 1 
10 43 13 5 zZ 2 1 


‘ihe primary tradeoffs that determine how rainout effects vary with 
inereased yiela are between (1) decreased scavenging probability arising from 
increased debris altitude and (2) increased debris concentration and 
debris-cloud size. The net result is a yield dependence of rainout effects 
that is qualitatively similar to the KAr dependence discussed in the previous 
section. 

The results in Figs. 42-45 and in Table V1 show that rainout effects have 
a broad maximun at yields of v2-5 kt (particularly at higher IWED) and a 
relatively rapid falloff below 1 kt (especially at lower IWBD). Rainout 
effects are also decreasing at 10 kt and will do so rapidly for higher yields 
(Sec. Vill-A-4). 

An important yield-dependent factor is evident in the maximum-threat- 


‘distance envelopes of Fig. 46. Primarily because of the larger fission-product 
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Fig. 46. 
Maximum dose-vs-distance envelopes for ten- 
shot salvos with RAF = 0.10. 


inventory, higher yield means longer threat lifetime at a given IBD and, thus, 
a larger maximum threat distance, r,.,. The higher wind speeds generally found 
at the higher altitudes associated with higher-yield weapons also increase the 
threat distance. lable VI1 and tig. 47 show estimated r,., vs yield for 
various IwED levels. Table Vfl also lists Vv, the average velocity of the 
debris-cloud centroid in the calculations, and the quantity Ymax/¥? which gives 
an approximate measure of threat lifetime. 
Cc. Wormalization to German Kain Lata 

ln this section, we use west German rain data to estimate probabilities 
that the rainout events depicted in figs. 34-47 will occur. We chose annual 
precipitation data here and estimated probabilities over a typical 1-yr period, 


although we could just as well have used seasonal data to estimate seasonal 


probabilities. 
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TABLe VIT 


MAXIMUM TBREAT DISTANCE, ro4,, AND AVERAGE DEBRIS-CLOUD CENTHGID 
VELOCITY, v, VS YIELD FOK 10-ShUiT SALVOS WITH RAF = 0.10 


IwbD = 10 rem IwBD = 100 rem 
YIELD v Pmax Pmax/V Pmax Pmax/V¥ 
(kt) km/h (km) _(h) 4km) a GD ee 
0.4 31 89 2.9 43 1.4 
1 35 130 3.9 66 1.9 
2 39 150 4g 78 2.0 
5 47 €25 4d 170 3.6 
10 55 255 4.6 220 4.0 
Ilwbb = 200 rem IwkBD = 450 rem IwbD = 1000 rem 
Prax Pmax/¥ Pmax Pmax/¥ Pax max /V¥ 
(xm) (h) (km) (h) (km) (hn) 
32 1.0 Ze 0.7 14 0.5 
Nb 1.3 28 0.6 17 0.5 
54 1,4 32 0.8 16 0.5 
120 2.6 T2 1.5 39 0.8 
210 3.6 170 3.1 115 2.1 


figure U-3 of Appendix G gives the climate-based conditional PbrF for the 
TEMPEST variable hAb,*® given that a precipitating storm system is in the area 
such thet he #0, where np is the 6-h precipitation accumulation. ‘he 


condition he #0, is equivalent to the condition HAr # 0, so we can write 


- 9 a 


*The large uncertainties in rig. G-3 stem directly from the hAlIW calculations 

of rig. Ged and will carry over pase into the final probability estimates. 
‘hese uncertainties are not just a result of poor statistics in the HAIN cal- 
culations; they are of more basic import. The statistical uncertainties in 

Fig. Gee derive poe from local inhomogeneities in the TemreSt-SCAWUEN 

rain field, on the spatial and time scales of the debris-cloud threats. Such 
large local variations occur in nature, and this fact injects a fundamental 
uncertainty element into any attempt to estimate, or even define, useful "“occur- 
rence probabilities" for rainout effects on the time and space scales pertinent 
to a specific tactical nuclear engagement. As we shall see, however, occurrence 
probabilities can still be boundea effectively. 
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Fig. 47. 
Maximum threat distance ise vs shot 


yield; ten-shot salvos with RAF = 0.10. 


Po os F(WAF > RAPO/RAF # 0) 
= probability that mAr > RAF 
= 1 = P( KAP < RAB /Re £0), 


or given that RAF # 0 


where P(HAF < KAY /ie # 0) is taken from Fig. G-3. Table Vil lists selected 


values of Pa. 


TABLE VIII 

RAF, Fy (vominal) Uncertainty bounds on Py 
0.015 0.97 0.92-1.00 

0.03 0.89 0.76-0.98 

0.05 0.75 0.56-0.94 

0.1 0.46 0.25-0.68 

0.2 0.16 0.06-0.26 

0.3 0.04 0.02-0.08 

0.4 <0.01 - 
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We need another type of climatic date, the probability that RAF # 0, to 
get overall probability of occurrence. Figure 48, a composite graph of West 
German data, shows the fraction of observation periods throughout the year that 
have observable precipitation versus length of the observation period. We used 
two different sources 9» 10 to get te, b=, and 24-h observation times. The 
plotted points reflect different stations available from each source and 
include stations in the Alps as well as in lowland and coastal sites. The 
"nominal" line appears adequate to represent the region of interest, with an 
urcertainty spread of perhaps +0.05 for times pertinent to our rainout 
scenarios (1-6 h). 

Wwe can use Fig. 48 to estimate probabilities of observable precipitation 
for different observation periods and thus make some allowance for the fact 
that the relevant "observation" period varies with yield and with IWED (Table 
VI1).* We define a) = PCR # 0) = probability that observable precipitation 
occurs during time T, and use Fig. 46 to estimate the following: 


*we should, in principle, also use different KAF Pbfs (tig. G-5) for different 
yields and different lwbbs. Such detail is presently not justified by the avail- 
able data, however, so we used the o-h PL throughout our work; this gives some- 
what conservative (nigh) ict d Peet estimates for our rainout events, which 
generally occur over times smaller than 6 h. 
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The overall probability that RAF # 0 during time T and that HAF > RAF is then 
.the product Py of P, (Table VIi1) and F.3 


Pp = P(HAF > RAFQ, Rp # 0) = Pyeo. 


Table 1X and fig. 49 show the resulting Pos Table X and Fig. 50 give the 
corresponding Pe range of uncertainty. 

Finally, we can identify the various scenarios (of Figs. 34-47) in terms 
of the probabilities in Tables IX and X (Figs. 49 and 50) by using Table VII, 
in conjunction with deposition-time information tabulated by ACRA, to estimate 
(conservative) maximum threat times of interest. Table XI lists these times. 


TABLi Ik 


ESTIMATED PHUBABILITY ThAT RAF # O DURING THREAT 
TIME T AND THAT RAF > RAF, 


Pr = P4Po 

BARS T= th 2h 3h Yh 54 6h 

0.015 0.16 0.23 0.27 0.31 0.33 0.35 
0.03 0.14 0.21 0.25 0.28 0.30 0.32 
0.05 0.12 0.16 0.21 0.24 0.26 0.27 
0.1 0.074 0.71 0.13 0.15 0.16 0.17 
0.2 0.026 0.038 0.045 0.051 0.054 0.058 
0.3 0.0064 0.0096 0.011 0.013 0.014 0.014 


0.4 <0.0016 <0.0024 <0.0028 <0.0032 <0.0034 <0.0036 
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F = PROBABILITY THAT RAF #0 OURING TIME T 


ANO THAT RAF 2 RAF, 
Fig. 49. 


Estimated occurrence probabilities for indicated 


storm and precipitation conditions. 
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Fe = PROBABILITY THAT RAF # O DURING TIME T 


AND THAT RAF 2 RAB, 
Fig. 50. 


Uncertainty bounds of occurrence probabilities 
for indicated storm and precipitation conditions. 
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TABLE X 


UNCERTAINTY BOUNDS FOR PROBABILITY THAT RAF # G 
DURING THREAT TIME T AND THAT RAF > RAF: 


Uncertainty Bounds for Pa 


RAP, T, = lh 2h 3h 4h Sh 6h 
0.015 0.12 ~ 0.19 0.19 - 0.27 0.22 - 0.32 0.26 ~- 0.36 0.27 - 0.39 0.29 - 0.41 
0.03 0.099 ~ 0.19 0.16 = 0.26 0.18 = 0.31 0.21 - 0.35 0.22 = 0.38 0.24 - 0.40 
0.05 0.073 - 0.18 0.12 ~ 0.25 0.13 - 0.30 0.16 - 0.34 0.16 - 0.37 0.17 - 0.39 
0.1 0.033 ~ 0.13 0.053 ~ 0.18 0.060 = 0.22 0.070 - C.24 0.073 + 0.27 0.078 - 0.28 
0.2 0.0078 ~ 0.049 0.013 - 0.070 0.014 - 0.083 0.017 - 0.094 0.017 - 0.10 0.019 - 0.11 
0.8 0.0026 = 0.015 0.0042 - 0.022 0.0048 ~ 0.026 0.0056 - 0.029 0.058 - 0.031 0.0062 - 0.033 
TAELE XI 


MAXIMUM THKEAT-T1mME BOUNUS T FOK 10-SHOT SALVOS Or FIGS. 34-47 


Yield WED > rem “200, 450 ~—*T000 
(kt) rem rem rem rem 
0.4 3 2 1 1 1 
1 4 2 2 1 1 
2 5 3 2 1 1 
5 5 4 3 2 2 
10 6 5 4 4 3 


to summarize, we use Table Al to estimate the appropriate "observation" 
time 1 for a standard 10-shot salvo with given shot yield and for a given IwWbD 
threshold. ‘This T then identifies the appropriate Py in Tables IX and x 
(rigs. 49 and 50) that applies to the given case. 

Despite the large uncertainty bounds on Pas Tables 1X-al and Figs. 49 and 
50 give probability information of direct value for assessing raincut threats. 
For example, the probability of any rainout threat from local rain fields, for 


our calculated scenarios, ranges from 10-20% at yields below 1 kt and IwbDs 
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2200 rem to a maximum of 30-40% for yields 22 kt and IWbDs >10 rem. Also, 
KAFs larger than those we calculated (0.45) are unlikely in the German climate. 
The worst-case rainout situations (highest occurrence probability and largest 
rainout effects) are those with RAF between 10.05 and 0.2 (Re : 0.1-1 mm/h), 
and this maximum is very broad. The rainout threat falls off rapidly at higher 
RAFS because of decreasing probability of occurrence and falls off rapidly at 
lower RAFs because of decreasing average rainout doses. These are fundamental 
2 features of rainout that will also apply to other scenarios. 
fable X11 shows RAF and Re for three different rain-field regimes. Table 
XLII lists some representative, specific, rainout-effects data from our 
’ ealeulated scenarios for the three regimes shown in Table XII, to illustrate 
different types of "casualty" information available from our results. 
xLII gives conservative estimates for the number of people N with IWBD > D, and 
associated occurrence probabilities for the 1-kt and 0.4-kt salvos of 
Figs, 34-41. 


D. Comparisons with Surface~burst Fallout and Prompt Collateral-Damage Effects. 


Results presented in the previous sections show that casualties could 
oceur from the postulated scenarios, and that although occurrence probabilities 
in the west German climate are not large they are not negligible. The actual 
number of people threatened with "damaging" doses in these engagements would 
depend on the level of radiation dose assumed to be damaging, but tens of 
thousands of people could be involved at levels of %100 rem and higher. 

In an actual engagement, tne number of people that would receive rainout 
doses at such levels would be much smaller than our results indicate because of 
shielding and postdeposition protective factors (Sec. V). Furthermore, even if 
we had made less conservative estimates for rainout "casualties," the 


assessment of their significance in actual operational contexts would be a 


TABLE XI1 
THREE RAlWePIELD REGIMES 


7 wualitative 6-h Average Rain-kate 
ae: Description | REF Ke (mm/h) 

Light <0.03 <0.4 

Medium 0.03 to 0.2 0.1 to 1 


eavy 20.2 >1 
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TABLE XIII 


NUMBER OF PEOPLE N WITH IWBD > Doe AND ASSOCIATED OCCURRENCE PROBABILITY P, FOR SCENARIOS OF FIGS. 34-41 


Light-Rain Field Medium-Rain Field Heavy-Rain Field 
Shot 
si D. (rem) P(N # 0) Ninax poe N a N ioe N 
1 200 0.25 + 0.03 20 000 0.04 <10 000 0.17 ~10 000 0.04 10 000 - 20 000 
1 450 0.17 + 0.03 9 000 0.03 <7 000 o.11 4 000 - 7 000 0.03 4 000 ~- 9 000 
1 1 000 <0.17 7 500 <0.03 <3 000 <0.11 1 000 - 2 500 <0.03 1 000 - 7 500 
0.4 200 0.17 + 0.03 14 000 0.03 <4 5000 Q.11 4 509 - 13 000 0.03 4 500 - 14 000 
0.4 450 <0.)7 11 000 <0.03 <3 500 <0.12 3 500 - 10 000 <0.03 3 000 - 11 000 
0.4 1 000 <0.17 2 000 <0.03 <400 <0.11 400 - 7 000 <0.03 200 - 9 000 


complex problem more appropriately handled by military scientists. We have 
thu: nade no attempt at such assessments here. We can provide some useful 
insight into the relative magnitudes of these rainout threats, however, by 
comparing them with corresponding fallout and prompt effects for the same 
seenarios. 

Fallout, prompt effects, and rainout effects for 10-shot salvos are 
compared in Figs. 51~54 for shot yields of 0.4, 1, 5, and 10 kt, respectively. 
These rainout scenarios are our standard airburst salvos (Sec. X), with RAFs of 
0.1 and 0.015; the fallout and prompt results are for the same scenarios with 
all shots fired as surface bursts.* io shielding or postdeposition factors are 
in any of these calculations. 

kainout effects compare least favorably for shot yields of v1 kt and less, 
where fallout effects are smallest, and at IwBDs <100 rem. In this regime, the 
worst rainout effects are higher than surface-burst fallout effects by a factor 
of 2 or less and are higher than the prompt effects** by a factor of 2 to 5. 

In real engagements, rainout "casualties" would compare more favorably than 
indicated because (1) rainout-specific, postdeposition, protective factors such 
as runoff, seepage, and avoidance possibilities would reduce real rainout doses 
as compared with real fallout and prompt doses, and (2) the probability of 
rainout (10-20%) is smaller than for fallout (100% for surface bursts) and for 


prompt effects (100% in all cases). 


*The surface-burst calculations were made by LASL staff member T. W. Dowler with 
the LACOMP fallout code. 


##Reeall that prompt effects will also be present in the rainout (airburst) sce- 
narios. 
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Fig. 51. 
Comparisons of rainout effects to surface-burst fallout and to 
prompt effects; ten 0.4-kt shots. 
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: Fig. 52. 
Comparisons of rainout effects to surface-burst fallout and to 
' prompt effects; ten 1-kt shots. 
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Comparisons of rainout effects to surface-burst fallout and to 
prompt effects; ten 10-kt shots 
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Fallout and prompt effects increase monotonically with shot yield, whereas 
rainout threats ievel off at shot yields of v2-5 kt and then begin to fall off 
(to zero) above V10 kt. Maximum potential rainout effects thus rapidly fall 
below fallout effects as shot yield increases above v1 kt, and they eventually 
fall below the level of prompt effects at shot yields somewhere above 10 kt. 
These features are seen in Figs. 51-54. Ina "light" rain field, potential 
rainout effects for our 10-shot salvos are comparable to or less than the level 


of prompt effects (for IwbD > +100 rem) for all shot yields. 


k. G&NERAL CONCLUSIGNS 

we have no "position™ on the commoniy asked question, "Is rainout 
important?" We believe the question is pointless unless the answer is an 
unequivocal "no," and such an answer is not possible. The only condition that 
would permit such an answer would occur if scavenging efficiencies for 
low-yield, airburst, nuclear debris were so low that only trivial amounts of 
debris could be brought down even in very heavy rain fields. This condition 
has not been proved as yet and may well never be. 

we must therefore assume that precipitation scavenging of low~yield 
airburst debris is possible, and this means that significant local rainout 
radiation fields are ana priori possibility if a rain field is in the area. 
because beth the number and intensity of such "hot spots" will increase 
monotonically as the number of shots increases, dangerouus rainout areas on the 
ground will be inevitable if the rain field persists and the number of 
low-yield shots is large enough. In this event, some of these not spots 
clearly could fall on nearby, downwind cities or friendly troops; and we would 
then havé a potential for rainout "casualties." 

The proper rainout "question" under these circumstances is then manifold 
and comolex: (1) where, when, under what conditions, to what extent, and with 
what probabilities ean rainout doses cccur; (2) what are the driving physical 
elements, how do they affect the answers to the above questions, and how can 
such insight be used to ameliorate undesirable collateral~-damage threats; and 
(3) what are the operational implications of the answers to all of the above 
questions in realistic military contexts, and does rainout raise important 


long-range issues that planners and decision makers should address. 
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We have dealt witn questions (1) and (2) in this report. In the rest of 
this section we make generalized observations that summarize some of what we 
nave learned about the magnitudes, sensitivities, and basic character of 
rainout threats. We also offer our present views on some long-range 
implications. 

The most striking result of our studies is tne consistent behavior of 
rainout effects as functions of the three driving independent variables: 
yield, average rain rate (RAF), and scavenging efficiency (measured by the 
parameter K). The dependence of rainout effects cn each of these variables is 
shown in rig. 55. These variables are independent in the sense that each of 
them produces the rainout behavior of Fig. 55 independently of the other two. 
Hainout threats are small, for example, if the yield, or the RAF, or the 
scavenging efficiency is small; furthermore, if any two of these variables have 
values that could produce substantial rainout, the dependence of rainout 


effects on the third variable is as shown in Fig. 55. 


Low TRANSIENT | MAXIMUM- THREAT 


THREAT — REGIME —j>———— | REGIME 7: 


RAINOUT 
EFFECTS 


A= YIELD, OR RAF, OR SCAVENGING 
PARAMETER Kk 


Fig. 55. 
Rainout effects vs yield, RAF, or scavenging efficiency. 
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The important features illustrated in Fig. 55 are as follows. 

1. There is a range of A between zero and a rough "threshold" (A,) where 
rainout effects are innocuous because A is too small to support significant 
rainout threats. 

2. In the transition regime between A, and Ap, rainout effects increase 
rapidly with increasing A as the intensity and extent of ground-deposited 
debris reaches significant levels. 

3. There is a broad, shallow maximum in effects versus A over a large 
range of A; effects are relatively insensitive to A over this range because of 
complex tradeoffs among factors that tend to cancel one another. 

4. A maximum threat level E,., can be identified within a factor of about 
2 in our calculations. 

5. There is a limiting value A, beyond which effects fall off rapidly to 
zero (for A = yield), or beyond which the values of A are no longer credible 
(for A = RA) or are of no interest (for A = kK). 

6. Rainout effects decrease slowly fron Emax as A approaches A3. 

Table XIV shows our current estimates of the three ranges in Fig. 55 for 
the three driving independent variables. 

The existence of the broad maxima in rainout effects versus yield, KAF, 
and k (item 3) is the most important of these features. This behavior gives us 


confidenee that the worst-case regimes have been adequately identified and that 
TALLE X1V 


BST1mATED VALUES OF A,, Ap, and 43 IN FIG. 55 
FUR The DRIVING VARIABLES YIELD, RAY, AND K 


Variable AS Ap A3 
Yield (kt) K0.1 0.4 to 1 “10 to 20 
RAE £0.01 0.05 to 0.1 0.4 
Kb “0.05 0.5 to 1 30 


®Rough guesses based on extrapolations of our results. 


scavenging coefficient A(n™!) is A = KKO° 8 | where H = 
instantaneous rain rate in mm/h. 
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no special sensitivities have been overlooked that cculd lead to unpleasant 
surprises. The breadth and snallowness of the KA and k threat regimes also 
mean tnat unpredictable fluctuations in local hAfs and uncertainties in 
scavenging efficiencies are relatively unimportant over broad ranges of values 
that are most important for rainout assessment. These are key points: they 
allow us to bound rainout threats and probabilities (as described earlier) 
usefully and with confidence. 

A second generalized observation based on our results is that "fine" 
structure, or "local" details*® in the time and space gradients of other 
physical variables are relatively unimportant, especially for multishot 
scenarios, provided certain integral conditions such as overall fission-product 
inventory and average rain rate remain essentially constant over the time and 
Space regions of interest. Examples of such variables include the initial 
debris-cloud source distribution; the size distributions, the internal- 
intensity distributions, and the individual ground speeds of rain cells; and 
the local distribution of ground-deposited debris as affected by local 
transport features. Variations in overall, net rainout effects as functions of 
such details generally appear to be less than the scatter produced by inherent 
Statistical fluctuations or by the heterogeneous population-weighting. 

A third important class of generalized rainout behavior is the direct 
correlation between dose level and distance and direction from the shots, as 
illustrated in rigs. 8, 9, 11, 26, 36-41, 46, and 47. Large numbers of people 
can be threatened with hign dose rates, but only if they are directly downwind 
from or relatively close to the engagement, and then only if a substantial rain 
field is in that same area at the same time as the debris clouds. These are 
conditions that should be predictable in a given conflict area for the 
pertinent few hours before and after an engagement. 

kinally, we emphasize again that shielding and other postdeposition 
factors would generally reduce the rainout lwBDs calculated in this study, 
perhaps tenfold or more. Furthermore, simple protective responses, 
particularly if preplanned, would provide additional large reductions in Iw8Ds 


for affected people; such responses might include avoidance of precipitation 


#Ey “fine and “local” we mean of tne same scale as debris or rain-cell sizes and 
threat lifetimes, and smaller. 


71 


accumulations, moving to protected places, or simply leaving the area of a 
rainout "hot spot." All these responses could be effective and would usually 
be feasible during the time required to accumulate large rainout doses. 
Specific operational considerations would also reduce rainout hazards; such 
considerations include city avoidance procedures, minimum total yield to do the 
job, and timely assessment of local wind and rain conditions. We believe that 
realistic inclusion of postdeposition factors, protective responses, and 
operational considerations would lead to the general conclusion that the 
scenarios calculated in our studies would not pose serious rainout problems. 

In summary, rainout is real; it must be considered an integral part of 
nuclear-weapon effects studies just as any other delayed weapon effect. Wwe 
believe debris from nuclear wars can be dangerous to nearby populations, even 
for airourst weapons; and to try to prove otherwise is futile unless unequiv- 
ocal new evidence should somehow materialize to show that scavenging 
efficiencies are very low. 

At the same time, we believe that the probability of producing large-scale 
rainout collateral-damage casualties is not large and that resultant effects, 
if they materialize, need not be prohibitive. Rainout does not rule out the 
effective use of small fission weapons in the airburst mode for defensive 
nuclear engagements. Rather, the airburst mode substitutes a low-probability 
threat of delayed effects for the certainty of delayed effects to be expected 
from other militarily equivalent modes; and rainout effects will generally be 
mere innocuous if they occur. Wwe expect that, if the character and magnitude 
of rainout threats are properly understood and are taken into account in 
realistic operational contexts, rainout damage may generally be within 
reasonable bounds. We velieve the results of earlier studies raise an 
unnecessarily severe specter in this regard. 

we emphasize, however, that we cannot be certain of this until more work 
is done (by the military community we hope) to define and evaluate conflict 
situations that provide the context within which rainout effects must be 
evaluated. Conflict scenarios, nuclear-weapon employment doctrine, and force 
deployment should be included. The next basic problem is to determine the 
effects of rainout on these elements. Until these effects are determined, the 


significance of rainout will remain in question. 
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APPENLIX A 
MESUSCALE VARITAbLLITY OF PRECIPITATLUl 


An important measure of uncertainty in precipitation scavenging is 
introduced through the variability of precipitation rates in time and space. 

In particular, when the variations are on a scale corresponding roughly tec the 
size of a contaminant cloud subject to scavenging, the morphology of the 
contaminant-precipitation system becomes a dominant consideration. Two 
limiting ideal cases that can help in scoping the problem are (1) a small 
contaminant cloud embedded in a large, homogeneous rain shield and (2) the 
inverse, a small rain shower within a large contaminant cloud. heal 
Situations, however, seldom present these simplified cases; and we must 
consider some concepts of interaction. ‘he simplest useful interaction to 
study is the morphological intersection of rain and debris clouds. The results 
of this type of study must eventually be modulated by the dynamic interactions 
of entrainment and mixing within the circulation systems associated with 
precipitation. 

For travel times of 1-10 h, dispersing clouds of atmospheric contaminant 
have horizontal areas of v1-1000 km@, Cursory examination of precipitation 
patches shows that they cover 10-1000 kn?, so the debris and precipitation 
occupy areas of equivalent size, and, in estimating interactions, we must 
account for their respective geometries. Atmospheric structure on the scale of 
current interest is best investigated by means of weather satellites and radar, 
Supplemented by in situ rain gauges to provide ground truth. Data collected by 
these observation tools have been examined in a number of studies. Although 
each study had different objectives, we can compile a composite summary of the 
mesoscale structure of precipitation in a "typical" occluding cyclone. With 
such a Signature available, we can use traditional climate data such as storm 
tracks and precipitation records for individual stations to calibrate a 
statistical simulation technique and deduce time and space patterns for 
precipitation in a selected geographic region. 

The picture that evolves is a multiscale structure with each scale 
exhibiting details on smaller scales, down to the basic unit of a precipitation 
cell. ‘the major large-scale structure nas been described by Boucher"’ and by 


Nagle and Serebreny. '¢ From their discussions, it appears that the occluded 
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midelatitude cyclone can be broken into four zones based on the character of 
precipitation, as follows. 

1. Continuous stratiform precipitation in a broad snield ahead of a 
developing wave (generally north to northeast). This shield covers a zone of 
a 104 km* in which the precipitation is distributed in a banded structure 
covering 10-25% of the zone. The precipitation is shallow and rates are low, 
leading to inefficient wet deposition of suspended material. this category is 
marked by the least amount of definitive observational data. sot all storms 
contain an area of extensive uniform precipitation. Summer and early-fall 
systems will have convective clouds and precipitation in tne northeast 
quadrant. 

2. tagged, patehy, stratiform precipitation in an area of 1000 km close 
to the storm center in the forward quadrant. This zone tends to represent 
patches of heavier precipitation superposed on the lighter background rates. 
Again, 10% of the zone is involved in active precipitation, with V>% at rates 
exceeding 6 mm/h. 

3. Frontal convective showers tnat are an integral part of mid-latitude 
eyclones. The frontal band during the active developing stage-of a storm may 
extend 1000 km to the south or southwest of the vortex center. The band tends 
to be long and narrow. with an east-west width of 50-100 km. The frontal band 
is not a solid band of precipitation but rather is represented by patches of 
precipitation covering vl0% of the region. tates exceed 6 mm/h in <1% of the 
zone. 

4, Postfrontal convective showers that tend tc be smaller, shallower, 
less intense, and more sparse than those within the frontal zone. These 
showers are more randomly distributed but still tend to crient themselves in 
bands. The postfrontal snower region covers a large area, but active 
precipitation in the zone probably covers well under 100 km< at any one time. 

tigure A-1 shows a composite of the models of Nagle and serebreny, !4 
Austin and houze, !5 and boucher, !! with characteristic radar echo signatures 
for each of tne four major regions from battan.!4 The figure indi*tes the 
inhomogeneity of precipitation echoes, but the absence of guantite rate 
information uncoubtedly results in an underestimate of the variations. 
Clearly, the size, deptn, and number of fundamental cellular units depend on 


location within the main system, 
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Fig. A-1 
Composite diagram of precipitation zones in a mid-latitude storm. 


The relationsnip between geograpnic location and type of precipitation 
encountered depends on local topography. upstream moisture sources, and the 
preferred tracks of synoptic storms. Although storms vary widely in 
development and motion, climatologists have been able to identify a few most 
likely storm tracks for a given region. ‘these tracks can be related to 
properties of atmospheric general circulation and orographic conditions. 

Figure 4-2 represents the two dominant storm tracks for western turope derived 
trom woronicz.© The northernmost track, which tends to be a summer case, would 
subject west Germany to frontal convective and postfrontal air~mass showers, \ 
The southern track would tend to oceur in winter and produce light, continuous 
precipitation. The reduced incidence of storms between the two dominant paths 
is due partly to the orographic barrier of the Alps. 

Gn the subsynoptic scale, many variables concerning precipitation 
structure are considered: rain area fraction (hAt), central peak rainfall 
intensity, intensity as a function of spatial position, cell size distribution, 


cell lifetime, cell shape, cell movement, orientation of cells' major axis. and 


a 
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Fig. A-2 
Major storm tracks over Western Europe. 


spatial variation of precipitation intensity. These quantities may vary within 
each synoptic zone above, but most of our calculations to date vary only the 
rain area fraction in each zone in the TEMPEST and SCAWUEN programs. Most of 
these quantities are given statistical distributions specified as cumulative 
probability distribution functions, 

The open literature has a number of fine papers summarizing various 
statistical results of radar precipitation patterns.''"@@ table A-1 presents 
selected parameters of precipitation patterns for each of the four major zones 
of an occluding cyclone. The table entries were derived from the literature 
references, and as many corroborating data sources as possible were used. ‘The 
entries reflect a measure of subjectivity in compromising among various 
reported results, but these compromises were generally easy to make because 
there was a remarkable consistency among the authors. The largest differences 
appeared in how the quantities were defined: different researchers defined 
Similar structure in either descriptive or statistical terms, or both. 

Gur simulation model is based on the parameters of ‘lable A-1, and the 
initial inputs were derived from the table entries. As the computational 
properties of the model evolved, however, we also recognized certain 
consistency requirements that changed the character of the parameterization. 
As the evolutionary nature of the model is important in understanding its 
present form, we will discuss some parts of Table A-1 and then present the 


- q@urrent parametric scheme. 


ty 
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TABLE A-I 


SUMMARY OF PROPERTIES OF MESOSCALE PRECIPITATION STRUCTURE 


(Numbers in square brackets denote references used in deriving table entries.) 


Property 


Area (km?) per Storm 
{11,12,13] 


Fraction of Storm Area 
with Precipitation 
20.25 mm/h 

[13,15] 


Maximum Cell Area 
(km?) 
{12,15,16,17,18,19] 


Altitude (km): 0.5 mm/h 
{13,14,15, 20] 2.5 m/h 
6.0 m/h 


Hourly Av Rate (mm/h) 
{15,16,21] 


Lifetime, T (h) 
{13,15,16] 


Shape 
[13,15,16] 


Orientation 
(with respect to W/D) 
{13] 


Movement 


16) 


Rate variation within 
cells 
{22] 


Continuous Patchy Frontal Postfrontal 
Stratiform Stratiform Convective Convective 
10° 10° 10° 10° 
0.10-0.2 G.10-0.2 0.10 0.02 
1 2 
f(A) = exp|——> (2nA ~ m) 
A 210 20 
m m m 
m= 2.5 to 3.2, 02 = 0.3 to 1.0 
3 4.5 5-10 5-6 
fe) 2 1-5 1-3 
fe) 1-2 0-3.5 0-1 
1 3 3 <0.1 
2/3 ee 
Scale-dependent: Tv 0.8 L (L = characteristic 
size). 


Height-dependent: 0.08 h per 2-km vertical extent. 
Ellipse, with aspect ratio @ given by the cumulative 
distribution function: 


8° 1.0 0.8 


P(B > B*) ie) 0.25 0.6 1.0 


Degrees, 8, measured clockwise, given by the cumula- 
tive distribution function: 


6° 90 100 110 


P(@< 6°) 0 0.15 0.70 0.95 


Approximately parallel to mid-tropospheric wind 


(Eq. A-3), 0.25 < y < 0.5 
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The area coverage is a function of the type of precipitation or the 
location in a storm. Table A-I shows that the continuous-stratiform region 
probably has the largest area covered by precipitation, although generally of 
low intensity. If we select a clock~hour rate criterion of 2.5 or 6 mm/h, the 
ragged patchy stratiform and the frontal-convective zones show the most coverage. 
zawadski 16 performed a detailed case study on a developing frontal wave that 
passed through Montreal on September 16, 1969. The gross statistics he calcu- 
lated from the radar returns can be combined to give a simple proportionality 


relation between the mean rainfall rate and its standard deviation: 
921.6 <k> , (A-1) 


wnere the angle brackets denote spatial average. ‘Zawadski used autocorrelation 
functions to deduce the mesoscale structure of the precipitation. The 
Lagrangian autocorrelation yields information on the characteristic lifetime of 
precipitation units. ‘Zawadski reports lifetimes of identifiable precipitation 
features to be about an hour. 

Tne two-dimensional autocorrelation yields a shape characteristic that is 
essentially elliptical with a ratio of major to minor axes of about 1.5 to 2.0, 
which is consistent with other studies. For the smaller scales (<10 km), the 
shape is nearly circular. 

The storm system moved with the individual cells until the final stage, 
when the storm appeared to be nearly stationary, with cells forming on the 
upwind side, moving through, then dissipating at the downwind edge. Although 
aawadski did not compare the movement with observed winds, the cells and system 
appeared qualitatively to move parallel tc the upper-level flow. 

Ine precipitation-cell size distribution derived from a composite of radar 
results of marshall and holtz,!> konrad et al.,25 and Austin and houze!> and 
from rain-gauge results of Fullerton and wilsoné4 appears to satisfy the 


log-normal relationship given in bq. (A-2): 
ee ll = 2 
£ (Ay) = oxy (tna, =)", (a2) 


where A, is the maximum cell area and m and 9, ere parameters of the 
ah as * * . : : : : ~oS ey - . 
distribution given by: median = e™, mode = e™-°li, ine chserved cell size 
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distributions demonstrate a great deal of consistency among the different 
geographical locations of tne observation sites. 

we recognized early in our work that processes such as precipitation 
scavenging that depend on the colocation of airborne debris and precipitation 
have a fundamental resolution requirement determined by the characteristic 
scales of variation of the two functions. in the present application, this 
resolution translates to a few kilometers and to a few tenths of hours. 
Conventional rain data do not provide this resolution, but, by the techniques 
discussed above, we were able to simulate precipitation with the required 
resolution. An acid test of such a model, of course, is its ability to match 
the observed conventional precipitation characteristics when they are 
appropriately integrated (see also Appendix F). The current parameters are the 
result of tuning efforts to produce precipitation that matches records for 
4emin, ji-n, and 6<-h rainfall observations in Germany. 

Several parameters described in Table 4-I that are factored into the model 
were determined to have a weak input toa either simulated rain-gauge records or 
scavenging-event probability. These parameters received less attention than 
those that showed greater sensitivity, such as tne cell size distribution, the 
distribution of maximum rain intensities, and the rain area fraction, which 
bear very directiy on precipitation integrals. On the other hand, the aspect 
ratio and the orientation of cells appeared to affect the integrated rainfall 
much less. Also, transitions from one zone of a starm to another had only 
Minor effect but complicatea the interpretability, so parametric studies were 
performed with 6- to 10~h traverses through a single zone. 

The results of KAér variation are shown in Fig. A-3. The RAF generally 
remains at 40% or less. Use of much larger values results in severe 
overestimates of the 1- to 6-h precipitation totals. The fhAF is basically a 
real physical parameter describing that fraction of a large area that 
precipitation is covering at a given time. The parameter is also related to 
cell size distribution and cellular-intensity cross section. 

in simulating a mesoscale precipitation structure, we must address several 
items of precipitation-rate information. Two main items are the distribution 
of peak rain rates for various cells and the form of rain-rate variation within 
the cell. bata for defining either parameter in detail are scarce, although 
existing data can certainly afford guidance in selecting logical functions and 


magnitudes. The fundamental precipitation units require rate data at very high 
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Fig. A-3. 


Effect of rain area fraction on distribution of 
rain rates. 


resolution to document their statistical properties. Two data sources are the 
l-min-resolution data of Jones and Sims¢4 and the work of Fullerton and 
wilson,®* which reflects resolution of a few seconds. These studies show 
precipitation rates that are very high compared with data that are smoothed 
over periods of an hour. Figure A-4 shows the distribution of peak rates 
derived by Wilson from radar and.in situ measurements. Kates exceeding 100 
mm/h are observed, whereas the median rate for cell cores is 25=30 mm/h. 

The decrease in precipitation rate from the core of a cell to its edge has 


been expressed in the form 


2 
B pes 
R(r,) = R, ad Beak ree ee i ee (A--3) 
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KONRAD et al. 


CASES (% BELOW INTENSITY) 


te) ; 60 160 
RAINFALL INTENSITY (mm/h ) 


Fig. A-4, 
Distribution of peak precipitation rates. 


There is no rigorous basis for selecting Eq. (A=3), and existing data 
offer little promise of imminent empirical documentation. The function does, 
however, offer a reasonable variation that has flexibility in the parameter Y 
to adjust "peakiness" and account for regions of very low background rates with 
narrow shafts of heavy rain (small y). When HAF, y, and cell sizes are all 
large, areas of low rain rate would be extensive and peaks of high rate would 
be few. As cell sizes go down, the number of peaks would increase. If yY were 
increased while the other two parameters stayed constant, the peaks would be 
broader, yielding more rain. As RAr goes down, the rain coverage drops from 
continuous to sparse. This method of describing rain morphology is thus quite 


general and includes such usual descriptions as “continuous light rain" and 


"scattered showers." 


APPENDIX B 
ACTIVITY NORMALIZATION AND INFINITE WHOLE-BODY DOSE (IWBD) BST”.MATION 


The SHASTA code estimates free-field IWBD from ground~deposited debris by 
means of the following assumptions: : 
1. Biological dose rate, d,, from deposited debris is proportional to the 
energy-decay rate of the debris; and this is, in turn, proportional to t°, 


where t = time from detonation and c is a "decay constant." 
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2. hainout-deposition times of interest, tas lie between a few tenths of an 


hour and a few tens of hours after detonation. 
53. The proportionality factor between dose rate and deposited-debris decay 
rate is a constant such that? 


d,; = dose rate from deposited debris at 1 h after detonation. 


2600 (rem/h) per (equivalent kt/mi2@) of deposited-debris intensity. 


Thus, 
d, = dose rate from ceposited debris at time t after detonation 


, 


= d4to ’ (B-1) 


and, for a sCAWUBLN-ShAS‘A simulation, 


= 2600¢°(rem/h)/(kt/mie) 


cr 
i] 


6734t° (rem/h)/(kt/km@) 


ee i 
= 6734Y JK t¢(rems/n) , (b-2) 


where 1 = fission-equivalent yield (kt), Liq = deposited-debris intensity 
(particles/km®), and ! = number of source particles in the simulation. 

tor gercss fission-product xctivity, the parameter c can be taken to be 
constant as a function of time, over the time span of assumption 2, with errors 


of less than 10 to 20%. Also, for such deposition times, IWbbD can be estimated 
within <10-30% by 


Bi opera 
IWBD = f d.dt , (B-3) 


4 
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provided x is a few nundred or larger and ty is in the range of assumption <. 


Bcuations &t-1 to b-3 tnen give 


l+c 
v (k -_1) 
ee Tee | teed 
d 
6734(ke*o - 1) | Thx 140 
2 . a 
ae 7 ty (rem/h) (B-4) 
In the calculations for this report we used c = -1.24 and k = 200 to get 
IWBD = 3d. t, = 2.02 x lo’y Hik tone 
i. Cine mM “da ‘ (B-5) 
APPbWD1d C 


‘bb NU&n wlub VATA bASt AL SPOUK TNTEKROLATIUN SCHEMES 


‘Ine built-in mean-wind data we use in SrlUun came from data tapes compiled 
by the wational Center for Atmospneric nesearch (WNCAH), boulder, Colorado. 
ihe NCAK compilations contain worldwide historical wind data translated to the 
fvational metecrolojical Venter (WMC) octagonal grid system.¢! tigure L=-1 shows 
the northern hemisphere of this system, wnich is a regularly spaced rectangular 
grid laid out on a polar stereograpnic projection of the eartn's surface. We 
later give tne transformations between HC (1,c) coordinates and geographic 
(latitude-longitude) coordinates. we have written codes* to read the WCAh 
tapes, process and print the data in various ways, tally and print statistical 
features of the data, and produce a SPOUk wild-data tape. 

The specific data base now in SPOOR covers a 42-point, western Europe 
subset (I = 33, 34, 35, 36, 37, 38; J = 23, 24, 25, 26, 27, 28, 29) of the WMC 


grid. Figure C-2 shows tnis subset. 


#This was the work of kh. Rich and k. Allenson, LASL. 
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Fig. c-l, 
NMC 47 by 51 grid. There are 1977 data 
points in the octagon. The pole point is 
I,J = 24, 26. 


The NCAR data are given in (1,J) space and consist of orthogonal wind 
components, WU in the I-direction and V in the J-direction. For each of the 42 
(1,J) grid points shown in rig. C~2, the specific data set extracted for SPOOR 
gives these components at 6 pressure altitudes (350, 706, 500, 400, 300, and 
250 mbar) for each 12-h interval from September 1965 to Lecember 1967.* 

Because SPOUK operates in a continuous phase space, we needed an 
interpolation scheme to supply WU and V values between the data-base points; and 
we developed a generalized code module for this purpose.”© When called by 
| SPOUR, this module generates biquadratic or bicubic least-square fits of the 
appropriate data-base winds versus I and J for any of the six altitudes and 


l2=-h time points in the data base. Linear interpolations in altitude and time 


then provide the required u; and Vj: 


*many of these data are not actually available. Wwe nave used simple linear 
extrapolations and interpolations to infer missing data-base entries. 
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Fig. C-2. 
NMC grid from SPOOR overlaid on map of Western Europe. 


illustrate, consider a trajectory step beginning at w 5 = (Ij, Jj Zs 
J 
ich defines bounding altitudes z, and z,,1 and times t; and ty,, in the 


se, so that 


£2; S$ 2iggd > 


make Separate least-square fits versus I and J over the 42~point grid 


C-2 for altitudes k and k + ' at times i and i+1. Linear 


interpolations in z and t between these surfaces then give u; = U(;) and Vj = 


Vw 5). 
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It is convenient to do these calculations in the (1,J) space, whereas the 
S'OOh tracking module operates more conveniently in the horizontal (x,y) space 
defined at each step by the tangent plane to the earth's surface at ® 5. A 
transformation between the (I,J) and (x,y) spaces is thus required. Wwe avoid 
some complications by using instead the transformation between (1,J) 


eoordinates and geographic coordinates. ©? As an example, for step j at Oa let 


the displacements in (x,y) space be (4x5, dy ;), where 
Ax., Ay.) = -4t, u.At * 
( x, y,) Ww, uy ) 

The corresponding displacements in (I,J) space are 
AI., 43.) = (Ay,/AE., Ax,/AE. 
( 4’ > (Ay, 4° y/ ? > 


where 


a, = 381/k. 
aE / 1" > 


. 1+ sin 60° 


Kk, 1 + sin %; 7 


5 = latitude of the point 
I.,J 
r, > ? 
ba” - Re 
sin ¢ 7) ; 
4a +R 


=, - 24)* + (I, - 26)- , 


4a = 973.712023868 . 


*in SPOOR, we arbitrarily take the positive x-direction along the positive 
Jeaxis and the positive y-direction along the positive i-axis. 
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Substituting the (1,J3)-to-geographical transformation for the 
(1,J)-to-(x,y) transformation introduces a small bias in the trajectories, but, 
at latitudes of usual interest to us, this bias is small for track lengtns up 


to several tinousand kilometers. 


APPLNDLX Lb 
Ra lN-MORPHOLOGY MODELING Tin THRE TEMPEST CODE 


‘tne rain-cell geometry models used in TbMPEST are based on the following 
assumptions (see also Appendix A): 

', The synoptic storm geometry does not change shape with time. 

2. All cells are elliptical and have a constant aspect ratio (ratio of 

minor to major axis) for a given cell as a function of time. 

3. Initial cell positions are uniformly distributed in area in each storm 

ZONE 

4, Linear cell-growth is parabolic in time. 

5+ The cell areas A,., at maximum size are log-normally distributed, 

ob. Lite-cycle times of nev rain cells are uniformly distributed on their 

individual lite-cycle time lines. 

These assumptions. the input Pls, and the logic flow described below 
completely aciine the TrmMrBS] methodology. Although some of these assumptions 
appear to be somewhat artificial. the parameters available in the formulation 
allow us to adjust overall cell structure and behavior to match statistically 
the characteristics of real precipitation systems. Furthermore, if we vary 
these parameters to study influences on final rainout effects, we find that the 
most important properties of the system are just those features that can be 
related empirically to real-world precipitation systems (for example, cell-size 
distributions, within-cell distribution parameters, and overall area 
fractions); these properties are exactly what is required of the model. 


To outline the model specifics, we refer to Fig. D=1 and define the 


following: 
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Fig. D-l. 
Rain-cell geometry. 


(X,Y) = TEMPEST coordinate system (defined by the storm-geometry input). 


(X_,Y_) = location of a rain-cell center. 


The 


where 
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¥ = (Vy9Vy) = rain-cell velocity (zone-dependent input). 
¢ = angle between the Y-axis and v. 
) 


= angle between the cell minor axis and v (obtained from an input 


PDF). 
=6@+o . 
= aspect ratio (obtained from an input PDF). 


= semimajor axis (see below). 


S oe BRB € 
i] 


= seuiminor axis = fa. 
equation for this cell is 


2 CXY+DX+EY+F=0 ,- 


+ BY 
ae cos" w+ sin’ WY > 
8° bin* vy + eos? yo, 
2 
2(1 - 8) sin pcos yp. , 
~2[x, (8° goa" U + eine W) + YG - 87) sin ¥ cos y] , 


~2[¥, (8° Sian y+ waar wv) + x, - 8°) sin » cos y] , 


(D-1) 


2 2,2 2 
F = x (8° oes ytsin py) +Y¥,(8 sin yt eee" wv) 


+ 2x YC - B”) sin wv cos J - Brae ‘ 


Both TEMPEST and SCAWUEN handle all rain cells in the form of Eq. (D-1). 
Because the parameters 6 and w are assumed constant throughout the lifetime of 


a given cell, the cell size is completely determined by the parameter a 


(semima jor axis). 
& cell's growth curve is given in terms of the life-cycle parameters in 


Fig. b-2 plus a correlation between a... and Tmay taken from the literature. 


anax = maximum value of semimajor axis during a cell's lifetime. 

Tmax = total lifetime of the cell. 

a _ > nes3 < 
Taax = 0.8 affe. (D-2) 


. 


When new cells are created (such as when TEMPEST sets up the starting cell 


structure), each new cell is assumed to start at a point in its life cycle 


selected uniformly over the interval (0, a aoe (assumption 6). To do this, we 


pick a 6T uniformly on (0, Tax) (Fig. D-2); the cell's size as a function of 


time is then 


=G +h + 172) 
a + + (D-3) 


MAX 


Fig. D-2. 
Cell growth curve. 
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where 


is 


I 


time measured from the cell's appearance (Fig. D2), 


© 


= NanaxdT (Tmax - 61) /Thax 


“Manax (Tmax - 26T)/Thax > 


“Nanay /Taax . 


we thereby completely determine the history of a cell. When we need a new 


eell, 
t. 


2. 
3. 


a value of A 4, is selected from an input log-normal distribution* 
(assumption 5): Eq. (D-2) gives 1, for the cell, 

a value of 6T is selected uniformly on [0, T,.,], and 

equations (D-3) and (D-1) then give the cell equations as functions of 


time. 


The overall TEMPEST logic flow is as follows: 
Inputs: Storm geometry, total leMPEST time for the simulation, cell veloc- 


1. 


De 


ities, rain-area fractions (RAts),** and PDFs for Away? 8: and @. 
Uniformly select an initial rain-cell location in a given zone of the 
storm system. 

Select values of A ..,, 8, and @ from the appropriate PDFs. 

Calculate coefficients for Eqs. (0-3) and (D~1). 

Kepeat steps 1 to 3 until the cumulative fractional rain-cell area in 
the given zone equals the input rain-area fraction. 


Repeat steps 1 to 4 until all zones are filled. 


This establishes the initial conditions for TEMPEST. 


6. 
ate 


Advance ‘all cells one time step in (X,Y), according to ¥. 

If, during any time step, a cell leaves. its zone or its area goes to 
zero kill it and start a new cell for that zone as in steps 1 to 3. 
Get new coefficients for all cells. j 

Repeat steps 6 to 9 from time step to time step until the total 


TEMPEST time is reached. 


This gives the desired TEMPEST storm and rain-cell data for SCAWUEN. 


#we use a routine given in Ref. 29. 
‘*ehAP is the same variable discussed in Appendixes A and G ard Section Ix. 
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APPEWNLIA E 
SPECIFIC MODBLS IN ‘THE SCAWLEN CUDB 


RALN-INTENSITY DISTRIBUTION WIThIN A PRECIP 1LTATION CrLL 
A. horizontal Distribution 

The horizontal (x,¥) distribution of rain intensity within a cell is 
assumed to be a time-invariant, two-dimensicnal Gaussian whose center intensity 
varies with time as the cell area and whose center value is Ko when the cell is 
at its maximum size. The parameter ho is picked from an input Por that can 


depend on the storm zone. Consider a specific cell at a given time step 


(Fig. t=1): 


(Xos¥Q) = location of the cell center (in SCAWUEN or SPUUn coordinates). 
(x,,y,) = location of a SPOCK particle within the cell. 

a = semimajor axis of the cell. 

b = semiminor axis of the cell = 8a. 

B = aspect ratio of the cell. 

w = angle between tne y-axis and the minor axis of the cell. 


n(%,d) = rain intensity at position (X,,¥g) in the cell. 


Fig. E-1. 
Precipitation-cell and SPOOR-particle 
geometry. 
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Then, 


~ 
] 


(x, - x) cos w + (y, - yi) sinw , 
d.= (y, - y,) cos w -. Ox, = 5) sin Ww , 


and the rain intensity (mm/h) at (x55¥,)> assumed to be proportional to cell 


‘area, is : 
; 2 2 2 
R(2,d) = [ 2 ). exp |-/£,4+4,}] , (E-1) 
a oO 2 “ : 
(Vmax 20 20. 
. a b 
- 
where 
a oe = cell semimajor axis at maximum cell size, 
o.* Yy2 ‘ 
C5 GPS 2 


The distribution parameters Ve and Yy are key input variables that can be used 
‘to match the ACRA system to rain-rate frequency and intensity data. 
B. Vertical Distribution* 

In our present rain~intensity module, we assume an effective rain intensity 
that depends on altiture but not on the other SPOOR or TEMPEST state variables. 
The Z~dependence is introduced through a simple multiplicative factor to account 
for the expected decrease in rain intensity with increasing altitude. The rain 
. intensity at time step j, corresponding to time t after the shot, at a location 
(x,»y,92,) within a rain cell, is given-by 


Ry (tgr¥geZ,) = Ry(ted) a(Z,) : (E-2) 


'.. where Rj (2.d) ds given by Eq. (E-1), and «(Z) is linearly interpolated from an 
input table.. 


*Also, see Appendix I. 


| i ; . 


4 


A more sophisticated model to acecunt for the variability of rain with 
altitude could readily be incorporated into SCAWULN. Such a model could be 
either deterministic or stochastic, and it could be a function of any of the 
ACKA state variables,* such as local winds, time, synoptic location in the storm 
system, and cell size. It cculd also use other types of input data such as 
PDrs for cell heights and vertical intensities. We have not inccrporated such 
a model because it is not clear that the sensitivity of overall rainout effects 
to such details warrants upgrading our system at this time. This question 


Probably deserves furtner investigation, however. 


11. SCAVENGING GF DEBRIS PARTICLES*® 

The scavenging rate for a SPOOR particle located at the phase-space 
position as (X55¥ 512554) in a local rain field with intensity given by 
Eq. (&-2) is 


aw 
aos: = - (R.) W , . (E-3) 
J s 


where Wy is the unscavenged debris weight represented by the SPOOR particle at 


> 
2s and Ay is a "scavenging coefficient." 
Typically, A’ depends on many factors and cannot be determined with high 


confidence. Sensitivity studies have established conservative values (highest 
rainout threat) for }, however; and its dependence on rain intensity can be 
included explicitly. The resulting expression that we use is 


k= 1.26 R9-8 (n-ty (E-4) 


where R is the rain intensity from Eq. (E-2). 
The debris weight scavenged by a rain cell from a contained SPOGK particle 


during a small time interval At is thus 


“rat 
Aw, = W(1 = eA) _ (E-5) 


> 
when a SPOOR particle at Qs is inside several overlapping rain cells, the total 


#Also, see Appendix H. 
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scavenged weight ‘during 4t is taken as a linear sum of the independently 


scavenged weights from each of the several cells: 


Sg > 
We = total weight scavenged during At froma given SPOOK particle at 25 
: ! 


i 
i 


ko ages 
~At > r. . 
J7l ii, A, At) ; (E--6) 


where iexl, 2, 3, ---, kK 1s a cell index for affected cells, the Ay are 
scavenging coefficients at a, for the individual cells, and k is the number of 
cells within which the SPOOh a eaeiéte is located at 25. 

After AW, has been determined for a SPOOR particle, a new particle with 
weight Aw, is created at 4, and the original SPUOK rarticle is continued to 
the next time step in its history with its weight reduced by AW, .# The set of 
newly created (scavenged) particles from all the SPOUn particles at each time 
step constitutes the input to the deposition and tallying modules. 


111. DEPOSITION OF SCAVENGD PAKTICLES#®# 
in SCAWJEN, each newly created scavenged particle is transported to the 
, ground as follows: 
Let v,, be the effective, vertical, water-transport velocity at altitude 2; 


then: 
= 11,2R9-°18 (msn), 
where 


K = rain-rate (mm/h) at 4 = a(2) Ky (Ry = ground rain rate). 


_ *When the SPOOR particle’ s weight decreases to a low value, the particle is 
subjected to a Russian Beutecce termination scheme, 39 at a substantial savings 
in computer’ time. 


r 


*#ALSO, see Appendix H. 


j 
& scavenged papticle located at vertical height ho (km) is then assumed to 


reach the ground ina time tp given by 


i! 


? Z Z 


6 
t.(h) = 92 = 0.09R ley —la__ * pe 
5 & fe) a(z)0-16 ny ) 


The ground rain rate is taken from kq. (E-1). During the time tr, each 
scavenged particle is transported houtzonte lly. according to a horizontal wind 
amen that varies linearly from a value VL, ) at a scavenging altitude Zo to 


¢, 


6.1 Vu ): at the ground; the horizontal velocity at any altitude « between the 


ground and z is thus given by 


a) = ¥(2 (> 1+ 0.9 Z| : (E~8) 


j 
The velocity ¥(Zo) is another input to SCAWUEN. Scavenged particles are 
tallied in the proper time intervals, as determined by te, and in a Lagrangian 
set of spatial bins on the ground see (See. VIII-A and below). 


> 
‘Typical v(2Z,) values range from perhaps 10 to 80 km/h. The horizontal 


displacements of scavenged debris particles before they reach the ground are 
typically <10-20 km; the more important rainout events, those with higher rain 
rates, imply horizontal displacements £10 km. Such displacements are worth 


including in ACnhA but do not warrant a more sophisticated treatment than that 


given here. 


1V, EFFECTS-TALLYING GRID* 

We assume rainout effects on the ground to be proportional to the 
intensity of locally deposited debris per unit area. This deposition intensity 
is determined by means of a Lagrangian tally-cell grid system erected at each 
time step in the calculation. This tally grid is formed by a set of annular 
rings centered on the ground-particle centroid, plus a set of K equally spaced 
angular divisions, as‘shown in Fig. k-2. The numbers of radial and angular 


#Also, see Appendix b and Sec. VILI-A. 
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Fig. E-2. 
Effects-tallying grid. 


bins are inputs to SCAWUEN and are fixed throughout tne calculation. Hadial 
dimensions of the tally bins vary with time and are based on statistical 
properties of the ground distribution of debris. 

To describe the radial binning for a given time step j, we define 


AW sm = weight deposited at step j from SPOOR particle m, 


(X5m9¥ jm) = ground location of AW jms 


M = total number of SPOUR particles. 


The centroid ‘of the AW5, at j is (Xjor¥ je)» where, if 


W = total weight deposited during step j 
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= z= 
W, 4 OMe ; (R-9) 
then 
M 
- > 
x. == AW. x, P 
je Wy foot ym jm 
M 
eee a) AW. Ys, (E-10) 
j m=] oh am 


fhe radial positions of the AW ins as measured from the centroid, are given by 


eae [ (*5m ~®_)° + (i = 7s he Geta 


and the first and second moments of P im about the centroid are, respectively, 


M 
eet D> Mira | : (E-12) 
and 
M 
3 are Ads 2 rad (B-13) 
= . md . . E= 
ty W, ae} jm" jm 3 


To get a radial binning structure at step j, SCAWUEN uses qs. (E-9)- 
(E-13) and the following recipe to set up the first N - 1 bin radii. 


radius of the nth annular ring from the center 


to} 
u 


1/2 


— 2 N -_ 
7 [22 eh —al n=1, 2, ..., N-l , (E-14) 
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where 


-r“/207 
d=e 4 > 
N = an input integer, 
ty = largest of the ae values. 


the remaining area, between ry_y and ry, is further subdivided into F 
adaitional annuli, each of equal area, where F is another input integer. 

This radial binning scheme is designed to give roughly equal numbers of 
deposited particles in each bin, a desirable feature from the standpoint of 
statistical-sampling efficiency; Eq. (E~14) is derived to give approximately 
equal areas in each radial bin by assuming a Gaussian distribution of deposited 


particles with variance One 4 
J 
The intensity of ground-deposited debris (particles/km@) in tally bin k at 


time step jis 


M 
A ; + 
I. = ao ; (E-15) 
J 4k 
where 
AW snk = weight deposited in bin k at time j by particle m, 
Ay = area of bin k at time 5° (xm). 


All tally-bin locations, dimensions, and debris intensities for all time 


steps are outputs to SHASTA, for final effects tallies. 
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APPENDIX F 
SPECIFIC MODELS IN THE RAIN CODE 


The xAIN code uses a standard TEMPEST output and executes a modified 


SCAWLEN simulation to tally rain rates and rain-frequency distributions at a 
fixed ground point. for a given time step At and a given rain cell passing 
over the fixed point during At, KAIN assumes the distribution of rain intensity 
on the ground to be given by the SCAWUEN intensity distribution for that cell 


over the given time interval (fq. E~1). RAIN then time-integrates this 


intensity over At, along the path traced out by the point across the (x,y) 


projections of the cell on the ground. This integration is performed for every 


céll passing over the point during each At in a TEMPEST-SCAWUEN calculation, 
and contributions from all cells are added together for each At to get total 


rain at the point during each At. These data can be binned and averaged over 


any other time interval (greater than At) to get appropriate rates and 

frequency distributions to compare with measurements. 
Figure i’-1 shows the geometry for a given integration. 

coordinates, consider a given time step j, and define the following for that 


Wie use SCAWUEN 


time step: 


(x,y) = SCAWUEN coordinates, 


(X-X0sY-¥Q) = translated SCAWUEN coordinates (as indicated in Fig. F-1), 


(X9,Y¥q) = location of the cell center at the given time step 


(from SCAWUEN), 


(Xp Vp) = location of the fixed ground point, 


V¥ = (vy,vy) = (constant) cell velocity, 
@ = semimajor axis of tne cell (from SCAwUEN), 
b = fa = semiminor axis of tne cell (from SCAwukn), 


P1,Po = intersections of the cell boundary and the extrapolated 
path of the point across the cell, 
(x',y') = coordinate system defined along the cell axes (fig. F-1), 


K(x',y') = rain intensity distribution for the cell (Eq. k=1). 
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PROJECTION OF 
RAIN CELL ON 
GROUND 


(%es Ye) 


PATH OF POINT 
ACROSS CELL 


Y-Ye 


Fig. F-l. 
RAIN integration geometry. 


The desired integration is then 


R = total rain during At at the fixed point from the given cell 


tz 


= h(x',y’) dt, 


ct 
~ 


where the integration is along the path of the point across the cell during At. 


The path of the point across the cell is a straight line with the equation 


v 
= es 4 aa 
y= ¥, + mx 12 °(3 Ye x) , 


or, in the (x',y') coordinate system, 
y' =A+ Bx' , 


where 
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x 
m= ra ; (F-4) 
y 4 
Ye = Yo = mx, ’ (F-5) 
y, 
yi at Oh ea (F-6) 


cos ¥ - m sin ¥ 


m+ tan ¥ 


B=jT-mtan¥Y- (F-7) 
The intersections Py and p. are given by the above equations, and 
a 
at ae _ ap \* - a? - b* (F-8) 
1,2 B* + 3? B* 3 ge Be + B ‘ 


we also neéd the location of the point p and the velocity, aye of p across the 


cell in (x',y') coordinates: 


= OS - sin Fm 
x x, cos Yp U4 (F-9) 
ve meee sin vw + v5 cos » , (F-10) 
V1 = “Vv, cos y+ Vy sinw , (F-11) 
Pp 
v., =-v. sint«-vcos;y . (F-12) 
Yp 


to effect the intcgration in Eq. (F-1), we note that the point travels 
across the cell at constant velocity (Vor, Vyt)s and we assume it to be at 
point p (see Fig. F-1) midway during the interval At. The times required to 


travel from P, to p and from p to pz are then 
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1 Vi, f 
x 
P 
x) - x! 
bt, ie Neeser ‘ (F-13) 
x! 
PB 


This establishes the integration limits for the interval At; namely, the 
integration is from t, = 6t, or At/2 (whichever is smaller) to to = Sto or 
At/2 (whichever is smaller). 

Because of the form of R(x',y'), ko. (F-1) cannot be integrated in closed 
form. 1t can be conveniently integrated, nowever, in terms of an approximation 


given by Hastings: 2! 


x 
; : 2 
4.[eeti-a hg (F-14) 
va JO 1 
a,x 
>i 
i=0 


where the a, are constants given in hef. 51. jhe result is 


i. 
p= ROOF jt  eegea sens eee | 1 er 1 ; (F-15) 


aly | c 4 4 
x (s ) 
a; 


Zn 


1=0 


where 
2 
F = AR . e) (F-16) 


, ; (F-17) 
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D= ‘ (F-18) 
D2... 22 
v8" + 1B 
2 
i Na 
B= Zi 5 : (F-19) 
v6" + Y\B 
Z, = ve (xs + 3) : (F~20) 
fy = ve (x3 + B) : (F-21) 


The parameters A, b, Yy, Yy» 8, and a are defined above and in Appendix E; the 


1 
when La and 4, have opposite Signs, and the negative sign applies when the 


signs are the same. 


APPENDIX G 
ESTIMATING The SPHCTRUM OF RAIN AREA FRACT1ONS 


The statistics of precipitation accumulation in ACRA depend on several 
parameters in the code system, including: 

1. Kain area fraction (RAF) 

2. Cell size distribution (oF, nm) 

3. Distribution of cellular peak precipitation intensity (R,) 

4, Cellular-intensity cross section (y). 

Synoptic meteorologists can estimate HAF in real time from weather-map 
analyses. This ability makes RAF a useful assessment parameter. RAF is also 
one of the most sensitive parameters for determining the statistics of 
precipitation accumulations. To minimize the degrees of freedom in modeling 


rainfall, we have tried to determine the other parameters for a given climate 


a; are the constants given in Ref. 31; the positive sign in Eq. (F=-15) applies 
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region and then let RAF be the variable tnat separates periods of heavy and 
light rain. This approach assigns the RAF a role of dummy, variable within the 
modeling context; we also feel that its interpretability as a physical quantity 
is what might be expected subjectively. 

Gur approach to tuning the model has been to use rainegauge records on 
various scales of time resolution, after scoping the range of parameter 
variation by using rainfall records as described in Appendix A. The model may 
be used to simulate the precipitation collected in a stationary rain gauge as a 
storm passes over it. This simulation has been made for several gauges 
separated spatially, with records of 0.1-, l= to 3-, and 6-n precipitation 
accumulations being simulated. The 0.1-h statistics were compared with 
,high-resolution (4-min) rainfall accumulations for two stations in Germany 
(Koblenz and Freiburg) so we could select reasonable values for all parameters 
except RAF. We then used 1- and 6-h accumulations to determine the spectrum of 
RAF values as follows: 

Figure G-1 shows the cumulative probability distribution of 6-n rain 
amounts. Figure G-2, produced by the model, estimates the corresponding 
accumulations as functions of RAF, by equating the accumulation statistic in 
Figs. G-1 and G-2, it becomes a straightforward matter to generate Fig. G-3 
(the cumulative density function of RAF, giver the existence of a rainy period) 
and the family of parameters describing cell size and intensity. Given the 
currently available data set, we feel fig. G-3 represents a reasonable spectrum 
of RAF values over tne middle 90%-frequency-of-occurrence range. 

we nave to address the rationale for equating the 6-h accumulations in 
Figs. G-1 and G-2 to construct Fig. G-3. Figure G-2 represents the rainfall at 
a number of sites for a given, statistically homogeneous region characterized 
by a specific RAF. Figure G-1 is also made up from data for a number of sites. 
The only criterion for counting an event, however, is that it must have rained 
at some time during the 6-h period; the properties of a precipitation region 
have no bearing on the count criterion. The selection of a time scale <6 h and 
the consistency between 1- and 6<h results suggest that, in general, equating 


the statistic with its model counterpart is justified. 
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Cumulative conditional probability distri- 
bution of RAF derived from climate data in 
Fig. G-1 and rain-case data in Fig. G-2. 
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APPENDIX H 
PRECIPITATION SCAVENGING 


J. Klett 


Debris particles in clouds become incorporated in the cloud liquid water 
or ice phases through nucleation and the action of a variety of transport 
mechanisms, such as convective Brownian diffusion, turbulent diffusion, 
diffusiophoresis, thermophoresis, electrostatic interactions, and relative 
gravitational sedimentation leading to hydrodynamic or "inertial" capture. 
Because only approximate quantitative descriptions exist for most of these 
processes, even when they occur in isolation and involve only idealized 
particles (for example, spheres or infinite cylinders with homogeneous physical 
properties), it is understandable why, as yet, no definitive assessment of the 


scavenging problem exists. : — 


ivevertheless, an important tentative conclusion may be drawn from the many 
idealized case studies that have been carried out in recent years, 32736 These 
studies imply that in-cloud scavenging by convective, precipitating systems may 
be regarded basically as a two-stage process. In the first stage, all the 
" scavenging mechanisms, and especially nucleation, brownian and turbulent 
diffusion, and diffusio- and thermophoresis, serve typically to transfer a 
major fraction of the debris particles (assumed here to be characterized by 
radiir x 1 um) to the cloud water, predominantly the fraction comprised of 
small droplets and/or ice crystals which together possess most of the total 
absorbing surface of the cloud water, in 20 min or less. In the second stage, 
this polluted cloud water is scavenged primarily through inertial capture by 
relatively large precipitating cloud particles such as raindrops, snowflakes, 
and graupel (soft hail) particles. This accretion process can cause a major 
fraction of the polluted water to fall out of the cloud in a time period 
similar to, though generally somewhat larger ‘han, that for stage one. Thus, 
net fractional depletions of debris-particle concentrations of the order of 
unity: may occur within the span of 1 hour, corresponding to an overall 
scavenging rate -d tn npdt of 1073 to 1074s"! (here n, and t respectively 
denote débris particle concentration and time). This overall theoretical 
" eonelusion is quite consistent with the observation, based on experiments in 
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the field, that in~cloud scavenging operates with approximately the same 
efficiency as the processes which convert water vapor to precipitation. 3? 

This description of the nature of in-cloud scavenging, in which the second 
stage is generally rate-limiting for the entire process, suggests that a 
reasonable upper bound for the effective scavenging rate is given by the 
accretion rate of small cloud particles by precipitation, with the former 
assumed to have already absorbed the debris aerosol through stage-~one 
processes. Clearly, such a formulation will also represent an upper bound for 
the case of below-cloud scavenging ("washout") as well, because the efficiency 
of the accretion process increases rapidly with the size of the collected 
particles. Experimental studies imply washout scavenging rates that are 
typically at least an order of magnitude smaller than corresponding in-cloud 
"rainout" rates, for the same precibitation flux.38539 We shall also ignore 
possible evaporation losses between eloud and ground. Although evaporation 
typically reduces the water flux at the ground to about half of that at cloud 
pase , 40 usually only a small fraction of this water loss is caused by drops 
that evaporate completely, allowing their captured debris particles to escape 
again into the environment before being "recycled" through further drop 
collection and breakup events. Finally, we shall “o:.sider only the case of a 
"warm" cloud with no ice phase, because terminal fall velocities and collision 
cross sections for most types of cloud ice particles are not well-known. 

Consider the collection of droplets of radius r' by raindrops of radius r 
>> r'. It is apparent that if there were no hydrodynamic deflection of 
approaching particles, the collection rate would be controlled simply by the 
relative velocity of approach and the geometric collision cross section. It is 
customary to describe the effect of hydrodynamic interaction on the collection 
rate by using the concept of the collision efficiency, EH = E(r',r), defined 
here to be the ratio of the actual collision cross section to the geometric 
eross section for a pair of interacting drops. Then if N(r,t) dr denotes the 
number of raindrops per unit volume on the interval (r,r + dr) at time t, the 


local depletion rate of the concentration n(r',t) of droplets may be expressed 


as 
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ones) = sate f K(r',x) N(r,t) dr , (H=1) 


min 


where Yinin denotes the minimum raindrop radius, and the "collection kernel" 


K(r',r) is given by 
K(r'yr) =™(r' #37)? E(r'yr)LV(r) - V(r')] . (H-2) 


) In this expression the qua.tity in square brackets is the difference between 
the terminal velocities of the drop and the droplet. We have also made the 
usual assumption that coalescence inevitably follows a collision event. 

Long41 has evaluated K(r',r) versus r for various radius ratios r'/r, 
using the collision efficiencies of Shafrir and Gal-Chen42 and Klett and 
vavis, 43 and the terminal velocities of' Gunn44 ana Beard and Pruppacher.45 An 
interesting result of Long's analysis is that K depends only weakly on the size 
ratio r'/r; for 10 < r < 500 um, K varies over about eight orders of magnitude 
for fixed r'/r, whereas for fixed r it changes by less than one crder of 
magnitude for 0.2 < r'/r < 0.9. ‘This implies that the scavenging rate 
coefficient, A(r',t), defined by the integral in Eq. (H-1), is only a weak 
function of r'. 

Let us now assume a steady state rain spectrum according to the empirical 


description of marshall and Palmer: 46 


N(r) = Noe rinem , (H-3): 


where No = 1.6% 10° m3 en! anda = 62 K70+21 em™', and where K is the rain 


rate in mm n7!, (Subsequent more detailed studies have shown that although 
pain ‘spectra often have exponential forms as in Eq. (H-3), the constant Ny is 
itself usually a function of Rf, and the functional dependence of a on K varies 
. somewhat from the above. 47 however, there appears to be no need for such 
refinements in the present application.) 

In view of the results of Long, 4! we may now combine kas. (h-1)-(h-3) to 


obtain the following simple estimate for the maximum possible scavenging 


“9 


- 


coefficient, which is independent of r' and t: 
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A ~w aN 5 f ve) ear . (H-4) 
v re} 3 


To obtain this expression we have made use of the fact that K is well 
represented by its form for r'/r << 1; E thus represents a characteristic 
collision efficiency for drop-droplet interactions, and is approximately unity. 
The replacement of r.i, by 0 in the lower limit of the integral is of little 
consequence, because of the factor r2¥(r) in the integrand. 

Also, from the observations of Beard and Pruppacher, 45 and considering the 
factors ré and e~°F in the integrand of kq. (h-4), we may represent the drop 


terminal velocity adeqvately by its variation over the range 50 < r < 500 um, 


namely, 


V(r) & Cr, 50¢ 4 ¢ 500 um , (H-5) 


where C is a function of air pressure and temperature. 
we may eliminate reference to the parameters Ny and C by introducing the 
definition of rain rate. When expressed as a mass flux of rain, the rain rate 


R' is 


alee ‘ 3 
R' = xr N(r) V(r) dr, (H~6) 
0 


3 


where a is the density of water, and we have assumed the absence of an 


appreciable updraft. If k' is given in cgs units (i.e., as g em72 s-'y, then 


the rain rate in mm n7' is kK = 3.6 x 104 R'/9 Therefore, from 


Eqs. (H-4)-(h-6) we have 


= 3 -or - 
E 0 e dr E a 
= = (z) ’ (H-7) 


A 

RY 478 4 

4.ax10° f[ 4 ~ar 4.8 x 10 
-jJgore dr 


which on substituting a = &2n79+21 becomes 
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1 


Ot. Regn 4 (H-8) 


a(s"*) 4 4.2 x 10.) ER 
It is clear from the derivation that the accuracy of this estimate should not 
vary with height 2, assuming the description N(r)/N, = exp (-8? R79*21r) holds 
for all z. 

A similar dependence of Aon k for the rain scavenging of large particles 
is evident from the early numerical calculations of Chamberlain,48 
Furthermore, excellent agreement between Eq. (H-&) and the results of a 
numerical computation by Crandall et al.34 of A(R,r') versus R for r' 2 10 cm 
can be achieved by choosing E = 0.63. Adopting this value, our final estimate 


for Ais 
A= 305 x 107R 7967) = 13k %n-) Rint (H-9) 


It 1s interesting to consider whether there are any conditions under which 
one might expect a significantly larger scavenging coefficient than Eq. (H=9). 
The one plausible situation that comes to mind is the thunderstorm environment, 
where strong turbulence and electrical forces might cause the rapid 
"self-collection" of drops of similar size, in addition to an enhanced 
aceretion rate of small drops by larger ones. however, the fact that rain 
spectra similar to the Marshall-Palmer distribution given by Eq. (h-3) are 
observed in thundershowers implies that a more rapid coagulation of cloud water 
would tend to be balanced by a larger drop-breakup rate, so that the overall 
Spectral shape would remain largely unchanged. 

Nevertheless, there appears to be a significant difference between the 
thunderstorm and nonthunderstorm cases, in that the stronger drop interactions 
in the former case imply the debris particles will be more uniformly 
distributed over the drop spectrum; i.e., the pollutant-to-cloud-water mixing 
ratio will be a weaker function of drop size. This in turn implies that the 
scavenging rate in a situation of strong drop interactions may be estimated 
Simply by the rate of total water depletion from the precipitating volume 
containing the debris. for example, if we suppose this volume is characterized 


by a vertical extent H and a horizontal-cross section A, and that its liquid 


410° 


water content is w (= water mass per unit cloud volume), then it contains a 
water mass of order Anw. If the rain rate is Ff", this amount of mass evidently 
will cross the bottom surface of the volume in a time T, where AR'T % AhW. °- 


hence, we estimate the scavenging rate to be Neg x qt! &% R'/bN, 
Note that the ratio K'/w is just the mass-weighted average velocity of the 


precipitation, <V>: 


m(r) V(r) N(r) dr 


0 (H-10) 


<Yy> = 


m(r) N(r) dr 
0 


where m(r) is the mass of a drop of radius r. This is the relevant effective 


transport velocity for the debris-laden rain. trom Eq. (h-3) we find 


2 4mN_9 cd = 8TN 9 
(ois 3) = anon ff ad — - 8.9 x 10° 7R9°84 (H=11) 
a 


which along with Eq. (H-10) leads to 
<v>(m o)) 23.2926 pan mm nt (H-12) 


This expression for the effective velocity of precipitation is consistent with 
data reported by Fullerton and wilson.49 {However, it disagrees considerably-- 


by about 60% for R = 1 mm h~!--with a result obtained by Kessler,29 who 


assumed, somewhat arbitrarily, that the velocity of precipitation is nearly the 


same as the terminal velocity of the median volume diameter particle [the drop 
with diameter D, such that half the water is on drops with L < Do, so that 


w/2 = [Pomc NC) aD].} 
with Eq. (h-12) we thus arrive at the following estimate of the scavenging 


rate under conditions of very rapid drop coagulation: 


0.16 
-1 
agls } Fy = Hinm . (H-13) 


lll 


For example, it H = 1 km and R = 1 mm h7', then Ay %*3x 1073 s-1, compared to 
A = 3.5 x 1074 s7’ from Eq. (H+9); i.e., Av. is about one order of magnitude 
larger than A. 

It may seem puzzling to note from Eqs. (h-¥) and (H=13) that A exceeds Ag 
if 70.6 2 104/h, This points up tne fact that for sufficiently large rain 
rates, the accretion of contaminated cloud water by rain occurs in a time 
period which is smaller than the time needed for fallout of the contaminated 
rain from the debris-containing cloud volume. For such rain rates, Eq. (h-13) 
provides a better (and smaller) estimate for tne scavenging coefficient. As an 


example, if kh = 1 km, then Eq. (H-9) becomes unrealistically large for k Pa 50 


moh ~!, 


' APPENDIX I 
RAIN@-LNTENSITY VERTICAL WEIGHTING FUNCTION 


The altitude dependence of precipitation intensity, designated by a(zZ) 
where Z@ is altitude, was derived from O'Reilly,?! who offers a family of curves 
for different ground-level intensities. when these are normalized to the 
surface intensity the various curves group together, so it is reasonable to 


select a Single function to represent 


a(Z) = —— , 


where R(Z) = rain intensity at altitude 4. This function is tabulated in Table 
I-1. -It shows a maximum at cloud base and decreases to zero at 19 km. The 
intricacies of intensity variation in the subcloud layer, and the cloud-to- 
cloud variations with altitude are not addressed here. 


112 


TAbLE I-1 
NOKMALIZED VERTICAL@INTENSITY WEIGHTING 
FUNCTION a(Z2) versus 24 


Altitude Zz 

__ (kn Z 
0 1.0 
0.5 1,0 
1.0 0.97 
1.5 0.90 
2.5 0.67 
3-5 0.39 
4,5 0.29 
5.5 0.16 
6.5 0.09 
7.5 0.03 
8.5 0.01 
9.5 0 
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